arXiv:1509.06606v2 [astro-ph.SR] 17 Mar 2016 


Astronomy & Astrophysics manuscript no. solan ©ESO 2016 

March 21, 2016 


Activity trends in young soiar-type stars ^ 

J. Lehtinen^, L. Jetsu^ T. Hackman^, P. Kajatkari^ and G.W. Henryk 

* Department of Physics, PO Box 64, 00014 University of Helsinki, 00014 Helsinki, Finland 
e-mail: jyri . j . lehtinenOhelsinki . fi 

^ Center of Excellence in Information Systems, Tennessee State University, 3500 John A. Merritt Blvd., Box 9501, Nashville, TN 
37209, USA 

Received date / Accepted date 


ABSTRACT 

Aims. We study a sample of 21 young and active solar-type stars with spectral types ranging from late F to mid K and characterize 
the behaviour of their activity. 

Methods. We apply the continuous period search (CPS) time series analysis method on Johnson B- and V-band photometry of the 
sample stars, collected over a period of 16 to 27 years. Using the CPS method, we estimate the surface differential rotation and 
determine the existence and behaviour of active longitudes and activity cycles on the stars. We supplement the time series results by 
calculating new log RJjj, = log F'hk/'^T^b srnission indices for the stars from high resolution spectroscopy. 

Results. The measurements of the photometric rotation period variations reveal a positive correlation between the relative differential 
rotation coefficient and the rotation period as A; oc but do not reveal any dependence of the differential rotation on the effective 
temperature of the stars. Secondary period searches reveal activity cycles in 18 of the stars and temporary or persistent active lon¬ 
gitudes in 11 of them. The activity cycles fall into specific activity branches when examined in the logProt/Pcyc vs. logRo^', where 
Ro“* = Iflr^, or log Prot/T’cyc vs. logRJjj, diagram. We find a new splif into sub-branches within this diagram, indicating multiple 
simultaneously present cycle modes. Active longitudes appear to be present only on the more active stars. There is a sharp break at 
approximately logRJjj, = -4.46 separating the less active stars with long-term axisymmetric spot distributions from the more active 
ones with non-axisymmetric configurations. In seven out of eleven of our stars with clearly detected long-term non-axisymmetric spot 
activity the estimated active longitude periods are significantly shorter than the mean photometric rotation periods. This systematic 
trend can be interpreted either as a sign of the active longitudes being sustained from a deeper level in the stellar interior than the 
individual spots or as azimuthal dynamo waves exhibiting prograde propagation. 
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1. Introduction 


Extended time series observations are among the most useful re¬ 
sources for studying the phenomena related to magnetic activity 
of stars. While the general activity level of a star can quickly 
be estimated from a single observation of its chromospheric or 
coronal emission level, most of the relevant phenomena are time 
dependent and require us to follow the stars over a longer period 
of time. Importantly, the variations of the observed activity in¬ 
dicators are often periodic in their nature ranging from the rota¬ 
tion signal to decadal cyclic activity variations. Thus, time series 
analysis is a crucial tool for studying stellar activity. 

Although rotation itself is not an activity phenomenon, it 
plays a crucial role in stellar dynamos. In particular, differen¬ 
tial rotation is an important parameter and determines whether 
the stellar dynamo more resembles an aQ or an dynamo 
(ICharbonneaull201^ . Quenching of differential rotation is pre- 


* Based on observations made as part of the automated astronomy 
program at Tennessee State University and with the Nordic Optical 
Telescope, operated on the island of La Palma jointly by Denmark, Fin¬ 
land, Iceland, Norway, and Sweden, in the Spanish Observatorio del 
Roque de los Muchachos of the Instituto de Astrofisica de Canarias. 

** The analysed photometry and numerical results of the anal¬ 
ysis are published electronically at the CDS and available 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsarc.u-strasbg.fr/viz-bin/qcat7J/A-l-A/yyy/Axxx 


dieted for rapidly rotating sta rs transitioning into dynamos 
(iKitchatinov & Rudigedll99^ . 

The magnitude of stellar surface differential rotation has 
been investigated with a number of different observational meth¬ 
ods. Most of them rely on the idea that dark starspots or bright 
chromospheric active regions forming on different latitudes will 
display different rotation periods and produ ce multiple varying 
periodic signals in the observed time series (lHalllll99Th . A mas- 
sive study based on t his p rinciple was recen t ly und ertaken by 
iReinhold et al.l (120131) and iReinhold & Gizoril (1201 5l) using Ke¬ 
pler data. Previous studies have revealed a positive correlation on 
the one hand between the relative differential rotation and the ro¬ 
tation period (iHalllllOOTh and on the other between the absolute 
differential rotation and the e ffective temperature (iBarnes et al.l 
l2005UCollier Cameronll2007l) . 

The rotational signal also reveals the longitudes or rota- 
tional phases of t he major spots or chromospheric active regions 
(iHall et al.ll2009l) . In many stars these concentrate on active lon¬ 
gitudes where most of the activity appears on one or two nar¬ 
row longitudes, whi ch can stay intact for decades (lJetsulll9^ 
iLehtinen et al.ll20rTh . Long-lived stable active longitudes can be 
interpreted as signs of non-axisymmetric dy namo modes present 
in the stars (lRadlerlll986HMoss et al.lll995l) . 

On longer time scales, time series observations using 
both photometry and chromospheric lin e emission have been 
used in the search of activity cycles (iBaliunas et al.l 119951: 
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Table 1. Basic observational properties of the sample stars, the comparison (Cmp) and check (Chk) stars used for the photometry, and the time 
span of the photometric record. 


Star 

HD/SAO 


B-V‘ 

d [pc]^ 

Sp. 

Cmp 

Chk 

photometry 

PW And 

HD 1405 

8.86 

0.95 

21.9 

K2V^ 

HD 1406 

HD 1439 

Aug 1988-Jun 2014 

EXCet 

HD 10008 

7.66 

0.80 

23.6 

G9V^ 

HD 10116 

HD 9139 

Dec 1998-Jan2014 

V774 Tau 

HD 26923 

6.32 

0.57 

21.2 

GOV^ 

HD 27497 

HD 26292 

Dec 1998-Feb 2014 

V834 Tau 

HD 29697 

8.09 

1.09 

13.5 

K4V-^ 

HD 284676 

HD 29169 

Sep 1993-Mar 2014 

V1386 0ri 

HD 41593 

6.76 

0.81 

15.5 

G9V^ 

HD 41304 

HD 42784 

Dec 1998-Mar 2014 

V352 CMa 

HD 43162 

6.37 

0.71 

16.7 

G6.5V^ 

HD 43879 

HD 43429 

Dec 1998-Feb 2014 

V377 Gem 

HD 63433 

6.90 

0.68 

21.8 

G5V^ 

HD 64467 

HD 63432 

Dec 1998-Apr2014 

V478 Hya 

HD 70573 

8.69 

0.62 

88.5 

G6V2 

HD 71136 

HD 70458 

Apr 1993-Apr2014 


HD 72760 

7.32 

0.79 

21.8 

KO-V-^ 

HD 71640 

HD 72660 

Dec 1998-Apr2014 

V401Hya 

HD 73350 

6.74 

0.66 

23.6 

G5V^ 

HD 73400 

HD 72412 

Dec 1998-Apr2014 

DX Feo 

HD 82443 

7.05 

0.78 

17.7 

KIV^ 

HD 83098 

HD 83821 

Feb 1992-May 2014 

FQHya 

HD 82558 

7.82 

0.93 

18.3 

KOV^ 

HD 82477 

HD 82428 

Nov 1987-Apr2014 

NQUMa 

HD 116956 

7.29 

0.80 

21.9 

G9V^ 

HD 114446 

HD 119992 

Dec 1998-Jun2014 

KUFib 

HD 128987 

7.24 

0.71 

23.6 

G8V^ 

HD 127170 

HD 126679 

Dec 1998-Jun2014 

HP Boo 

HD 130948 

5.86 

0.58 

17.9 

F9IV-V^ 

HD 128402 

HD 127739 

Dec 1998-Jun2014 

V379 Ser 

HD 135599 

6.92 

0.83 

15.6 

KOV'^ 

HD 136118 

HD 137006 

Dec 1998-Jun2014 

V382 Ser 

HD 141272 

7.44 

0.80 

21.3 

G9V^ 

HD 141103 

HD 139137 

Dec 1998-Jun2014 

V889 Her 

HD 171488 

7.39 

0.62 

37.2 

G2V2 

HD 171286 

HD 170829 

Apr1994-Jun 2014 

MV Dra 

HD 180161 

7.04 

0.80 

20.0 

G8V2 

HD 182735 

HD 177249 

Feb 1999-Jun2014 

V453 And 

HD 220182 

7.36 

0.80 

21.9 

G9V^ 

HD 219224 

HD 221661 

Dec 1998-Jun2014 

V383 Fac 

SAO 51891 

8.57 

0.85 

23.7 

KIV^ 

HD 212072 

HD 212712 

May 1994-Jun 1996 
Oct2011-Jun 2014 


References. iMontes et al.l ( l2001hh iGrav et al.l ( l2003l iGrav et all(l200il iGrav et alJ Jloolh iGaidos et al.l ( l2000h 


iMessina & GuinanI 12002) : lOlah et al]|2009t) . This field of study 
is biased because the existing observational records are still 
quite short compared to the longest activity cycles that have 
been found. Moreover, the cyclic activity variations are typically 
quasiperiodic rather than stationary processes, which further de¬ 
creases the efficiency of the period search methods. Neverthe¬ 
less, it has been possible to relate the estimated cycle periods to 
other stellar parameters. This has revealed a sequence of activity 
branches, suggesti ng different dynamo modes being excited on 
different stars (e.g. lSaar & Brandenburg! 19991) . 

In this study we perform an analysis of 21 young nearby 
solar-type stars using ground-based monitoring photometry 
gathered over a period of 16 to 27 years. Most of the st ars have 
been selected from the sample studied bv lGaido^(ll998h . but the 
observing programme also includes six additional stars with no¬ 
tably high levels of activity. 

The stars in our sample are listed in Table[T]along with their 
basic observational properties. The Johnson V band magnitudes 
and B - V colours ar e taken from the Hipparcos and Tycho cat¬ 
alogues (IESAI IT997h . and represent the mean magnitudes and 
colour indices from the observations over the full Hipparcos mis¬ 
sion. The distances d are likewise derived from the Hipparcos 
and Tycho trigonometric parallaxes. The spectral types reported 
for the stars are taken from a number of sources given individu¬ 
ally for each star. 

The colour magnitude diagram of the stars is shown in Fig.[T] 
based on the Hipparcos and Tycho data and assuming negligible 
interstellar extinction. The plot a lso indicate s the zero age main 
sequence (ZAMS) according to ICoxI (l2000l) . It shows that the 
stars all lie on the main sequence although there is an apparent 
offset to larger absolute magnitudes for the coolest K-type stars. 


2. Observations 

2.1. Photometry 

Our study is based on photometry obtained with the T3 0.4 m 
Automatic Photoelectric Telescope (APT) at the Fairborn Ob¬ 
servatory in Arizona, which has been monitoring our programme 
stars since late 1987. In this paper we include all of the standard 
Johnson B- and V-band photometry collected with the telescope 
up to June 2014. The time span of the observations is given in 
Tabled and ranges from 16 to 27 years depending on the star. 
The only exception is a 15 year gap in the observations of SAO 
51891. 

The observations from the stars consist of differential pho¬ 
tometry, where the variable target stars (Var) are compared to 
constant, usually F-type, comparison stars (Cmp) as the differ¬ 
ence Var - Cmp. In addition the constancy of the comparison 
stars is simultaneously monitored by observing separate constant 
check stars (Chk) as the difference Chk - Cmp. 

We estimate the typical error of the target star photometry 
to be between 0.003 and 0.004 mag nitudes based on monitoring 
constant stars with the same setup (lHenrvlll995h . Errors of the 
check star observations can be expected to be somewhat larger 
since fewer individual integrations are used to determine their 
values. For a brief descri ption of the ope r ation of the APT and 
reduction of the data, see iFekel & Henr\l (l2005h and references 
therein. 

2.2. Spectroscopy 

To determine the chromospheric activity level of the sample stars 
we have observed their visible spectra with the high resolution 
fibre-fed echelle spectrograph FIES at the Nordic Optical Tele- 
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Fig. 1. Colour magnitude diagram of our sample stars. The plot symbols denote spectral types: yello w triangle fo r F-type, light orange circle for 
G-type, and dark orange square for K-type stars. The ZAMS is indicated by the grey line according to lCod ll20QOh . 


scope (iTelting et al.ll2014l) . This instrument provides full spec¬ 
tral coverage within the wavelength interval 3640 - 7360 A 
in 79 overlapping orders. The spectroscopic observations were 
obtained in 2012 and 2014. We performed the observations in 
the high resolution mode, which gives a spectral resolution of 
R - 67000. The observations were reduced using the FIEStool 
pipeline. 


3. Time series anaiysis of the photometry 

For the time series analysis of our photometry we used the con- 
tinuous period search method (hereafter CPS) formulated by 
iLehtinen et al.l (1201 ill . The method models the light curve data 
with non-linear harmonic fits, 

K 

%ti) ^ M + ^[Bk cosiklnfti) + Ck sin(k27rfti)], (1) 

k=l 

to short datasets selected from the full time series data by apply¬ 
ing a sliding window of predetermined length. The order K of 
the harmonic fits is adaptive and determined from the data using 
the Bayesian information criterion. In this study we allowed fits 
of the orders K = 0 (i.e. a constant brightness model), K - I, 
and K - 2. For each dataset we calculated the mean magnitude 
M, the full light curve amplitude A, the period P = /“', and the 
epochs of the primary and secondary light curve minima fmin.i 
and fmin, 2 - Naturally fmin.i can only exist for K - 2 models and for 
K - 0 models only M is defined. Fach of the fits was checked for 
reliability by comparing the distribution of the fit residuals and 
the error distributions of the model parameters against Gaussian 
distributions. If any of them was found to be significantly non- 
Gaussian, the whole dataset was labelled as unreliable. To find 


the initial search range for the light curve p eriods, we first ap- 
plied the three stage period analysis method (l.letsu & Peldl 19991 
hereafter TSPA) with a wide frequency range before proceeding 
with the CPS analysis. 

The lengths of the individual datasets are defined by a max¬ 
imum time span ATmax- If a dataset starts with a data point at to, 
all following data points within [fo, fo + ATmax] are included in 
the dataset. Because of the sliding window approach, most of the 
datasets have data points in common with the adjacent datasets, 
which will introduce inherent correlation in the model parame¬ 
ters. To overcome this correlation, we chose a set of independent 
non-overlapping datasets as the basis of our further statistical 
analysis. We also discarded all datasets having njata < 12 data 
points because their fits are likely to have low quality. 

As opposed to pr e vious st udies using the CPS method 


dFehtinen et al.l 


_ 1201 IL 120 121 iHackman et all l201lL 120131: 

iKaiatkari et aTTl20l4 l2015r we chose to use both the B and 
V bands in this study to get more data points for the periodic 
fits and to increase their precision. As can be seen in Fig. |2] for 
HD 171488, there is a linear correlation between coeval photom¬ 
etry in the two photometric bands. This applies generally to the 
photometry of all our stars. The linear relation means that the 
two bands contain essentially the same information. It is thus 
reasonable to simply scale one band on top of the other and use 
the resulting combined time series as a single photometric band. 
We used the V-band as the basis and found the empirical scaling 
relation B - ciV - cq between the two sets of differential pho¬ 
tometry separately for each star. The proportionality coefficients 
Cl are listed in Table |3 In each case we found ci > 1 with the 
mean value (ci) = 1.29, meaning that the B-band light curves 
always have larger amplitudes than the V-band curves. This be¬ 
haviour is consistent with modulation caused by spots that are 
cooler than the surTounding photosphere. 
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Fig. 2. B-band vs. V-band differential photometry for HD 171488 with 
a fit showing the linear correlation between the two colour bands. 


The use of combined B- and V-band photometry has in¬ 
creased the number of reliable period detections for all stars. For 
some of the lowest amplitude stars in our sample we were able 
to detect periodicity in 10-20% more datasets in the combined 
data than when using the V-band only. With the largest amplitude 
stars the use of combined B- and V-band data removed all = 0 
models, allowing period detection throughout the data. Using the 
combined photometry also increased the overall ability of CPS 
to find reliable model fits for the datasets. As a median, the num¬ 
ber of reliable fits was 8% higher for the combined data than for 
the V-band only. 

The crucial parameter to be set for the CPS analysis is the 
upper limit of the dataset length ATmax- To find the most reason¬ 
able value, we performed the CPS analysis for our sample stars 
with values ATmax = 20 d, 30 d, and 45 d. We found that the dis¬ 
persion of the period estimates, AP^ (see Sect. 14.11 Eqs.|2]and 
O, from the individual datasets was relatively small for the two 
longer values, but for ATmax = 20 d there was more scatter. This 
indicates that at Ar^ax < 30 d the dataset length and the number 
of data points contained within these datasets have become inad¬ 
equate and it is no longer possible to produce repeatable period 
estimates for many of our sample stars. If we intend to use the es¬ 
timated period fluctuations as a measure of differential rotation 
for the stars, we should use ATmax ^ 30 d. With these dataset 
lengths the estimated period fluctuations appear to be approach¬ 
ing a lower limit likely governed by physical processes on the 
observed stars. 

We can also define a time scale of change for the light curves, 
Tq, and see how this changes with varying ATmax- This is the 
time from the start of a dataset duri ng which a model fit c an ade¬ 
quately describe the following data dUehtinen et al.l20ril) . As ex¬ 
pected, we found that as ATmax increased so did Tq- This shows 
that as the datasets get longer the light curve of the observed star 
has more time to evolve and smear out finer details. When details 
are lost, the light curve fits tend to get simpler and consequently 
better describe future data. 

We settled on using ATmax = 30 d for all of the sample stars. 
This is a good compromise between getting reasonably precise 
period estimates and restricting the light curve evolution within 
the datasets. It is also similar to the dataset length s we have used 
previously i n other studies using the CPS (e.g. iLehtinen et al.l 
1201 lLl2012t[Kaiatkari et all2015ll . 



1995 2000 2005 2010 

year 


Fig. 3. Differential B- and V-band photometry from HD 171488 {top 
two panels) and the constant check star HD 170829 used for it {bottom 
two panels). The same magnitude scale is used for all the panels. 


Fig. |3] shows as an example what the differential photometry 
looks like for HD 171488. This is a star for which we see a clear 
periodic signal. As can be seen from the raw photometry, the en¬ 
velope of the light curve of HD 171488 indicating the light curve 
amplitude due to rotational modulation is clearly much wider 
than the scatter seen in the photometry of the constant check star 
HD 170829. Consequently, the CPS was able to find good peri¬ 
odic fits both for the V-band alone and for the combined B- and 
V-bands. There were, however, slight improvements associated 
with using the combined bands; the number of datasets with re¬ 
liable fits increased from 517 to 585 and the number of those 
where a period could be detected from 513 to 585. 

When using the combined bands, the variations of the esti¬ 
mated periods showed an expected drop with increasing ATmax- 
For HD 171488 using ATmax = 20 d, 30 d, and 45 d we 
found APw = 0.048 d, 0.012 d, and 0.009 d respectively. From 
ATmax = 20 d to ATmax = 30 d there is a large drop in APw indi¬ 
cating increased precision in period determination but increasing 
the dataset length to ATmax = 45 d left the level of precision prac¬ 
tically unchanged. This is an encouraging result if we intend to 
measure physical period variations on the stars. 

For the same three values of ATmax, we found the estimated 
time scale of change of the light curve to be Tc - 22.7 d, 28.3 d, 
and 36.9 d. There is a clear increase in Tc when moving from 
ATmax = 20 d to ATmax = 45 d, indicating the loss of finer detail 
in the light curve when using longer datasets. This shows that 
using longer datasets will lead to some level of information loss 
but also that the calculated time scale of change Tq seems to 
depend mostly on the chosen analysis parameters and not the 
observed stars themselves. 

Another star, HD 43162, with a lower light curve amplitude 
shows clearly the advantages of using the combined photometric 
ban ds. In our previous stu dy of the V-band photometry of this 
star (iKaiatkari et al.ll2015h we were able to detect periodicity in 
192 of the 370 datasets with reliable model fits. Here, with two 
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2000 2002 2004 2006 2008 2010 2012 2014 


year 

Fig. 4. Comparison of period estimates for HD 43162 from plain V- 
band photometry (top) and combined B- and V-band photometry (bot¬ 
tom). The P estimates from the independent datasets are denoted by the 
diamonds with error bars and from the rest of the datasets as points. The 
grey shading indicates the time span of all reliable model fits. 


Table 2. Properties of the light curve fits. The proportionality coefficient 
Cl between the B- and V-band photometry, the numbers of independent 
reliable K = 0 and K > 0 order fits, and the mean V-band amplitude of 
the K > 0 order fits in the CPS results. 


Star 

HD 1405 
HD 10008 
HD 26923 
HD 29697 
HD 41593 
HD 43162 
HD 63433 
HD 70573 
HD 72760 
HD 73350 
HD 82443 
HD 82558 
HD 116956 
HD 128987 
HD 130948 
HD 135599 
HD 141272 
HD 171488 
HD 180161 
HD 220182 
SAG 51891 


Cl n[c=o 


1.18 

0 

1.52 

30 

1.28 

48 

1.09 

0 

1.25 

1 

1.63 

14 

1.27 

7 

1.30 

3 

1.23 

23 

1.37 

21 

1.22 

1 

1.20 

0 

1.27 

0 

1.31 

36 

1.39 

8 

1.25 

24 

1.25 

0 

1.23 

0 

1.38 

15 

1.22 

3 

1.18 

0 


nK>o 

<A) 

77 

0.102 

16 

0.010 

15 

0.013 

52 

0.096 

70 

0.026 

48 

0.019 

68 

0.022 

100 

0.043 

52 

0.014 

40 

0.017 

130 

0.063 

117 

0.093 

81 

0.033 

35 

0.014 

66 

0.018 

39 

0.018 

63 

0.031 

82 

0.067 

47 

0.013 

55 

0.032 

15 

0.086 


more years of data, we found periodicity in 229 of the 407 reli¬ 
able datasets from the V-band photometry. This increased to 326 
out of a total of 425 reliable datasets when using the combined 
bands. Graphically this is visible in Fig. |4] where both the num¬ 
ber of period estimates and the total number of reliable model 
fits, indicated by the grey shaded area, increased when combin¬ 
ing the two photometric bands. Clearly the availability of more 
data points improved our capability to detect periodicity from 
the very noisy lowest amplitude light curves. 


4. Data analysis and results 

Here we describe the methodology used for analysing the photo¬ 
metric and spectroscopic observations. We summarize the basic 
characteristics of the CPS fits in Table by giving for each star 
the number of independent reliable datasets with non-periodic 
K = 0 and periodic K > 0 order fits as well as the mean V- 


band model amplitudes from the K > 0 fits. In Table [3 we then 
present the main results from the period analysis. Results from 
the spectroscopy are presented in Table |4] 


4.1. Rotation 


The photometric periods derived from the independent datasets 
for each star show various levels of repeatability from star to 
star. These period variations may be the result of sparse data , 
low amplitude, or low S/N observations (iLehtinen et al.ll20Tlh . 
On the other hand, they may be signs of surface differential rota¬ 
tion or active region growth and decay occurring simultaneously 
at various latitudes and longitudes on the star. For example, a 
large active region may be forming in one location while an¬ 
other is decaying at a second location well separated in longi¬ 
tude. The resulting shifting phases of light curve minima may 
result in measured photometric periods that do not correspond 
to the true rotation period on any latitude. We have been careful 
to minimize these problems by suitably defining the independent 
datasets as described in Section[3]above, in particular by limiting 
th e duration o f any dataset to 30 d. Therefore, we will assume, as 
in iHalll (Il991l) . that the observed range of periods for a particular 
star is a measure of the differential rotation with stellar latitude. 

We characterize the mean rotation periods of our stars by the 
weighted mean 


r, _ 

‘w — ^ 

LWi 

and the weighted standard deviation 


( 2 ) 


APw = 


ZWi(Pi - Pw)2 




(3) 


of the independent period estimates, where the weights are the 
inverse square errors of the individual periods, w, = cr'^. The 
mean rotation periods of the stars are listed in Table |3] 
where their errors are given as standard errors of the mean, 
cTp w = , calculated from the weighted standard devi¬ 

ation and the total number nK>o of independent datasets with a 
periodic model. As the differential rotation e stimate we use the 
relative +3cr range of the period fluctuations (I.Tetsulll9^ 


6APw 

Fw 


(4) 


As mentioned above, there are certain caveats to directly in¬ 
terpreting the Z value as a measure of differential rotation. If we 
can be certain that the entire range of periods is due to differential 
rotation, the Z value should approximate the relative difference 
between the fastest and slowest rotating spot areas on the star, ap¬ 
proximately corresponding to the 99% interval of the estimated 
period values. Moreover, if these extreme period values corre¬ 
spond to the rotation periods at the equator and the poles of the 
star, there would be an unambiguous connection between Z and 
the differential rotation coefficient k = AQ/Ggq- In practice we 
know from solar observations and Doppler imaging of other stars 
that this is generally not the case, and the observed spot distribu¬ 
tion is confined to a narrower latitude band. The observed range 
of rotation periods from the spots is thus expected to be smaller 
than the full range, thus predicting Z < ktme- On the other hand, 
observational errors, sparse data coverage, and the need to use 
short datasets to get local period estimates will cause additional 
uncertainty in the period values and cause somewhat increased 
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Table 3. Results relating to the weighted mean photometric rotation periods P„, the active longitude periods P^i, and the cycle periods Pcyc- 


Star 

Pw [d]“ 

Z 

Pal [d]“ 

Qk 

^cyc [yr] 

cycle type* 

FAP 

FAP grade* 

HD 1405 

1.75622(62) 

0.019 

1.752212(71) 

3.82- 10“" 

8.0 

MAM_M+ 

3.0- 

10-* 

good 

HD 10008 

6.78(13) 

0.47 

6.80 

0.81 

10.9 

MM_M+ 

1.1 ■ 

10-3 

fair, long 

HD 26923 

11.08(20) 

0.42 

10.66 

0.99 

7.0 

MM_M+ 

3.8- 

10-3 

poor 

HD 29697 

3.9651(59) 

0.065 

3.9433(28) 

9.24 ■ 10“^ 

7.3 

A 

3.8- 

10-3 

poor 

HD 41593 

8.135(30) 

0.19 

8.0417(23) 

3.45 ■ lO-*’ 

3.3 

M_ 

8.2- 

10-3 

fair 

HD 43162 

7.168(38) 

0.22 

7.1323(10) 

1.88- 10-3 

8.1 

A 

5.1 ■ 

10-3 

poor 

HD 63433 

6.462(40) 

0.31 

6.46414(98) 

4.20 ■ 10-3 

2.7 

MM+ 

1.2- 

10-3 

fair 






8.0 

AM_ 

2.5- 

10-3 

poor, long 

HD 70573 

3.3143(34) 

0.061 

3.29824(36) 

2.69 ■ 10-3 

6.9 

MM_M+ 

5.1 ■ 

10-3 

good 

HD 72760 

9.57(11) 

0.50 

9.68 

0.70 






HD 73350 

12.14(13) 

0.40 

12.59 

0.58 

3.5 

MM_M+ 

3.6- 

10-3 

fair 

HD 82443 

5.4244(43) 

0.054 

5.41471(22) 

2.50 ■ 10^^ 

4.1 

MAM_M+ 

9.2- 

10-3 

fair 






20.0 

M 

3.1 ■ 

10-3 

good, long 

HD 82558 

1.60435(42) 

0.017 

1.6037330(68) 

1.44- 10-3 

14.5-18.0 

MAM_M+ 

< 10-'3 

excellent, long 

HD 116956 

7.860(16) 

0.11 

7.84203(72) 

2.06 ■ 10-'* 

2.9 

MM_ 

4.0- 

10-3 

fair 






14.7 

A 

5.5- 

10-3 

good, long 

HD 128987 

9.80(12) 

0.44 

9.63 

0.25 

5.4 

MM_M+ 

4.7- 

10-3 

fair 

HD 130948 

7.849(26) 

0.16 

7.89 

0.85 

3.9 

MM+ 

3.3- 

10-3 

poor 

HD 135599 

5.529(68) 

0.46 

5.13 

0.99 

14.6 

MM+ 

2.9- 

10-** 

good, long 

HD 141272 

13.843(84) 

0.29 

13.80 

0.27 

6.4 

M 

9.1 ■ 

10-3 

poor 

HD 171488 

1.3454(13) 

0.054 

1.336923(41) 

3.98 ■ 10-3 

9.5 

MM_M+ 

2.1 ■ 

10-'* 

good 

HD 180161 

9.91(11) 

0.46 

9.41 

0.92 






HD 220182 

7.678(23) 

0.13 

7.62002(21) 

2.37 ■ 10-4 

13.7 

MM+ 

2.9- 

10-4 

fair, long 

SAO 51891 

2.4179(41) 

0.039 

2.40 

0.98 







Notes. Numbers in parentheses show the error in the last digits of the period estimates. *** Cycle type; "M" cycle found from M results, "A" cycle 
found from A results, "M_" cycle found from M - A/2 results, "M+" cycle found from M + A/2 results FAP grade: "excellent" PAP < 10“®, 
"good" 10“® < FAP < 10“^, "fair" 10“^ < FAP < 10“^, "poor" 10“^ < FAP < 10“*, "long" cycle longer than half the length of the photometric 
record 


values of Z. In general, the best scaling of the period variations 
into a differential rotation estimate, minimizing both over- and 
underestimation of k, is not an obvious choice. We recommend 
that our Z values be primarily used to determine the functional 
relation of the differential rotation to other astrophysical quanti¬ 
ties or as an indicator of whether the differential rotation of a star 
is weak or strong. We note that the exact numerical factor cho¬ 
sen for the calculation of Z does not affect the results presented 
in Sect. 16. II 

Because of the above mentioned uncertainties, we report the 
raw Z values in Table [3] There is a tendency for the large Z 
values to be found in stars with poorer observational coverage 
and fewer available independent period estimates. However, at 
the same time, the star with the poorest observational coverage, 
SAO 51891, has one of the smallest values of Z. In Sect. 16. II we 
further investigate the effect of the spurious period fluctuations 
on the reliability of differential rotation estimation from the Z 
values. 

4.2. Active longitudes 

The longitudinal distribution of major spot areas is tracked by 
the light curve minimum epochs fmin. Using the rotation period 
P, we can transform these into minimum phases ^min = (fmin - 
fo) /P mod 1, where to is an arbitrary epoch chosen to dehne the 
zero phase. Each dark spot on a star will contribute a depression 
in the observed brightness and for equally strong spots with a 
phase separation larger than A(/> ^ 0.33 these can b e observed as 
separate light curve minima dUehtinen et al.ll2()TTh . 


If there are active longitudes present on a spotted star, we ex¬ 
pect to observe a long-te rm phase co herence of fmin with an ac¬ 
tive longitude period P^i (|Jetsulll9^. This period can be found 
for example by using the Kuiper test (iKuipeJ 19601) with a range 
of folding period values. We applied t he unweighted Kuiper test 
as formulated by iJetsu & Peltl (1 19961) on the primary and sec¬ 
ondary light curve minimum epochs fmin.i and fmin ,2 from the 
CPS results. For each detected Pai we calculated the error es¬ 
timate using bootstrap. 

The active longitude periods and their critical level values 
2 k for the Kuiper test, i.e. the p-values against the null hypoth¬ 
esis of uniform distribution of fmin, are both listed in Table [3] 
We found active longitudes for nine of our sample stars with 
Qk < 0.01 and two more with 2 k < 0.1. These stars show a 
wide range of different behaviours from the very weak and only 
occasionally appearing phase coherence on HD 63433 to the sta¬ 
ble and persistent active longitudes on HD 116956. In the case 
of stars for which no Pal could be found with 2k < 0.1, we dis¬ 
play the best candidate period from the Kuiper periodogram in 
the vicinity of their P„. We have not calculated error estimates 
for these periods since they are very uncertain. 

As is visible from the light curve minimum phases plotted for 
the individual stars in Figs. BI-ITOl using their Pai, the estimated 
active longitude periods are generally only average values over 
the full observation records. Most stars with active longitudes 
show either some jumps in the minimum phases or more gradual 
phase migration. The best dehned migration patterns are seen 
on HD 70573, HD 82443 and HD 220182. These are discussed 
further under the individual stars in Sect.|5] 
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4.3. Activity cycles 

From even a cursory inspection of our plots below, it is evident 
that the active stars show long term variability that can be linked 
to variations in the activity level and o ften appears to follow 
cyclic patterns. LikelRodonoeLalJ (|2Q( ^, we applied the period 
detection method of lHorne & Baliunasl(ll986l) (hereafter the HB 
method) to four parameters from the CPS results; the mean mag¬ 
nitude M, the light curve amplitude A, and the combined values 
M - A f2 and M + A f2. These parameters describe respectively 
the axisymmetric and non-axisymmetric parts of the spot distri¬ 
bution and the minimum and maximum spotted area on the star. 
We excluded SAO 51891 from this cycle search since its pho¬ 
tometric record is dominated by the long gap of 15 years in the 
observations and the remaining data covers a time span that is 
too short to warrant reliable cycle detection. 

We grade these cycles by the false alarm probabilities {FAP) 
given by the period d etection method following the scheme of 
iBaliunas et alJ (Il995h : “excellent” for FAP < 10 “good” for 
10-^ < FAP < 10 ^ “fair” for lO^^ < FAP < 10 ^ and 
“poor” for 10“^ < FAP < 10 '. Our secondary grade label 
“long” differs from their scheme. It is given to cycle periods that 
have lengths larger than half the length of the whole photometric 
record and are thus less certain than the shorter cycle periods. 
We do not state the error estimates for the cycle periods given 
by the HB method since the true uncertainties of the cycles are 
dominated by their non-periodic behaviour and are impossible to 
estimate from the currently available data. 

We found cyclic behaviour in 18 out of 20 stars and evidence 
of two separate cycles in three of the stars. However, many of 
these cycle periods either have a large FAP or are of the same 
order as the length of the photometric record. In all except five 
cases the cycle periods are found in the M data. Typically the 
same periods are also present in the M - A/2 or M -H A/2 data, 
often with a somewhat larger FAP value. Overall, the A data 
shows much less evidence of cyclic behaviour than the other pa¬ 
rameters M, M - A/2, and M H- A/2. A cycle in A can be found 
only for seven stars. For three of them this cycle is unique to the 
A data. 

In the case of HD 82558 there is a strong signature of a long 
cycle in all four parameters. They fall into a range from 14.5 yr 
to 18.0 yr, suggesting that they correspond to the same physical 
cycle, but they do not give a clear indication of what the exact 
length of this cycle is. In this case we have to conclude that the 
cycle of HD 82558 is too long to allow more than an estimate of 
the period. 


4.4. Chromospheric indices 

We quantified the chromospheric activ ity of our sample sta rs by 
measuring the Mount Wilson S-index dVaughan et al.l[T97^ 


5=0-- 


H + K 
R + V 


Table 4. Chromospheric emission indices log/?Jjj, and activity classes 
from this work and the chromospheric emission indices 
from the literature sources used for calibration. 


Star 

logRjjK 

activity 

class'* 

log^HK,cahb 

HD 1405 

-4.217 

A 


HD 10008 

-4.480 

A 

-4.530^ 

HD 26923 

-4.618 

MA 

-4.555^ 

-4.521^ 

HD 29697 

-4.036 

VA 

-4.066^ 

HD 41593 

-4.427 

A 

-4.456^ 

HD 43162 

-4.425 

A 

-4.480^ 

HD 63433 

-4.452 

A 

-4.424^ 

HD 70573 

-4.488 

A 

O 

1 

HD 72760 

-4.609 

MA 

-4.454^ 




-4.57-^ 

HD 73350 

-4.700 

MA 

-4.657^ 

HD 82443 

-4.286 

A 

-4.264^ 

-4.234^ 




-4.30^ 

HD 82558 

-4.079 

VA 

-4.03-^ 

HD 116956 

-4.366 

A 

-4.447^ 

HD 128987 

-4.505 

MA 

-4.439^ 

HD 130948 

-4.533 

MA 


HD 135599 

-4.462 

A 

-4.663^ 

HD 141272 

-4.566 

MA 

-4.452^ 

HD 171488 

-4.175 

VA 

-4.21-^ 

HD 180161 

-4.541 

MA 


HD 220182 

-4.388 

A 

-4.503^ 

SAO 51891 

-4.327 

A 

-4.05-^ 


Notes. Activity classes: “VA” very active stars with logfJJjj, > -4.20, 
“A” active stars with -4.20 > logl?(jj, > -4.50, “MA” moderately active 
stars with -4.50 > logi?([j. > -4.75. 

References. iGrav et ^ i l2003ll iGrav et ^ il2006h IWhite et"^ 
l (2007h 


since they provide uniform measurements for most of our sam¬ 
ple stars. To improve t he calibration at th e strong emission end, 
we also used values of lWhite et al.l (l2007h . For normalized FIES 
spectra we obtained a - 19.76. 

The S-indices are transformed into fractional fluxes /?hk = 
FHK/mT^g in the line cores with respect to the bl ack-body lu¬ 
mino sity of the stars using the conversion formula (IMiddelkoonl 
h982h 


lO^'^Ccf^ 


(6) 


(5) 


Rhk = 1.34 
where the colour dependent conversion factor 
logCrf = 0.25{B -V)^- 1.33(B - Vf + 0.43(8 -V) + 0.24 (7) 


and transforming it into the fractional emission flux log /?' at 


the Can H&K lines. The process involves measuring the emis¬ 
sion flux through two triangular bands FI and K with FWHM of 
1.09 A centred at the two line cores and normalizing it to the flux 
at two flat continuum bands V and R with a full width of 20 A and 
centred around 3901 A and 4001 A. The normalizing constant a 
is needed to adjust the measured values to the original Mount 
Wilson HKP-1 and HKP-2 spectromete rs. We chose to calibrate 
our measurements against the values of iGrav et al.l (120031 1200^ = /?hk - ^phot- 


is applicable to main-sequence stars with 0.3 < B - V < 1.6 
(lRuttenlll984l) . This value still contains a photospheric contribu¬ 
tion which is described as 


logRpho, = -4.898 -H 1.918(B - V)^ - 2.893(B - Vf 


( 8 ) 


for stars with B-V > 0.44 (iNoves et al.ll984ll . This is subtracted 
from the /?hk value to get the corrected value 


(9) 
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Table 6. Age estimates adopted for the kinematic groups. 


Group 

abbr. 

age [Myr] 

AB Dor Moving Group 

AB 

50-120^ 

Her-Lyr Association 

HLA 

120-280^ 

Hyades Supercluster 

Hya 

600^ 

IC 2391 Supercluster 

IC 

35-55^ 

Local Association 

LA 

20-150^ 

Pleiades Moving Group 

Pie 

100^ 

UMa Moving Group 

UMa 

300^ 

Tuc-Hor Association 

THA 

20^0^ 

TW Hya Association 

TWA 

4-12^ 


Referen ce s. ^^*lMal o et al.lil2 bl.3h ^^^lEvans et alJi2012h *^^ lMontes et all 
( l2001bll <^> IZuckerman & Sond (l2004li 


The final logarithmic log values are given in Table 0] 
together with values from the sources used for calibrating the 
measured S-indices. All of our stars are either “active” or “ very 
active” according to the classification of lHenrv et al.l (Il996h . i.e. 
their log/Jj^j, > -4.75. To help differentiate our stars accord¬ 
ing to their chromospher i c activ ity level, we divide the “ac¬ 
tive” class of lHenrv et al.l d 19961) into two subclasses (“active” 
and “moderately active”) at its mean log/?^j^. The resulting ac¬ 
tivity classes are “very active” for log/?^j^ > -4.20, “active” 
for -4.20 > log/?^j^ > -4.50 and “moderately active” for 
-4.50 > log/?^j^ > -4.75. 

5. Individual stars 

In this section we discuss the noteworthy results for each of the 
sample stars individually and compare our results with past re¬ 
search. The main characterizing parameters from the literature 
are the rotation period estimate, the chromospheric Can H&K 
and the coronal X-ray emission indices log/Jj^j, and log/?x, 
the reported age estimates, and the identifications of kinematic 
group membership. These parameters are listed in Table 5 along 
with references to their sources. Ages adopted for the various 
kinematic groups are given in Table |6] This table also lists the 
abbreviations used for the groups. 

The results from the CPS analysis are presented graphically 
for each star in Figs. ISUTOl These plots show the M, A, and P 
results from the independent datasets for each star, as well as the 
primary and secondary minimum epochs fmin folded into min¬ 
imum phases (p^i^ with the best periods found by the Kuiper 
method. In the cases where the folding period does not have a 
significant Kuiper statistic and cannot be identified as an active 
longitude period, we have coloured the fourth panel in the plots 
grey. 

If not mentioned otherwise, the stars lack observed stellar 
companions. For some of the stars distant companions have been 
found and we note the basic characteristics of each of them. In 
all of the cases the observed companions are on wide orbits and 
the primary components in our sample are thus effectively single 
stars. 

5.1. PWAnd-HD 1405 

HD 1405 (PW And) is an “active” rapidly rotating (P^ - 
1.7562 d) K2V star. It is among the most active stars in our 
sample, with the emission index log/?^j. = -4.217. It is known 
to have strong rotationally modulated chromospheric emission 


(iMontes et alJ 12001 at iLonez-Santiago et al.l l2006t IZhang et"^ 
20151), and the emission in dex log/?[jj. = -3.85 reported by 
Lopez-Santiago et"^ (1201 Ol) puts it firmly into the “very active” 
class. The star is very young; isochrone fitting and identifica¬ 
tion as an AB or LA member place its age between 20 Myr 
and 150 Myr. The h igh lithium abundan c e is consistent with 
a Plei ades-type age dMontes et al.ll200T^ . IStrassmeier & Ric^ 
(l2006l) presented a Doppler imaging temperature map indicating 
that the spot activity is concentrated at latitudes below -1-40°. 

We found both a well-defined 8.0 yr activity cycle and active 
longitudes on HD 1405. As expected, the activity variations are 
not stationary and although the cycle has repeated itself multiple 
times during the observational record, none of the consecutive 
cycles have been identical. 

The active longitudes show more stability and the main activ¬ 
ity area has stayed for the most of the time at a single rotational 
phase in the reference frame of Pal = 1.75221 d. On some occa¬ 
sions, most notably before 1990 and between 2000 and 2002, the 
main activity area was located for a while on the opposite side 
of the star with a roughly A<p = 0.5 phase separation from the 
typical active longitude phase. The occasional phase jumping of 
th e main active area is analogous to the "flip-flop" events found 
hv Lletsu et al.l (1 19931) for FK Com, where the light curve minima 
alternated between two active longitudes on opposite sides of the 
star. During 1989, we saw one flip-flop event underway; during 
the course of the second half of the year the light curve minima 
quickly but gradually migrated from one active longitude to the 
other. 

5.2. EXCet-HD 10008 

HD 10008 (EX Cet, HIP 7576) is an “active” G9V star with age 
estimates ranging from below 10 Myr up to 440 Myr. Most of 
the age estimates, including the identification as a HLA or LA 
mem ber, are consistent with an ag e above 100 Myr. On the other 
hand, iNakaiima & Morin^ (1201 2l) identified the star as a TWA 
member implying an age below 10 Myr. The age discrepancy 
makes the TWA membership seem unlikely. 

HD 10008 has a low amplitude light curve. We were able 
to detect periodicity in only 16 of the 46 independent datasets. 
The mean V-band light curve amplitude in the periodic K > 0 
order fits is O'l'OlO. We found evidence of a 10.9 yr activity cycle 
in the photometry. There is an apparent anticorrelation between 
M and A such that the light curve amplitude is highest when 
the mean magnitude is at its faintest. This indicates that the spot 
activity becomes increasingly non-axisymmetric on this star as 
the activity level rises. Nevertheless, the HB method was unable 
to find a cycle period from the A results. 

5.3. V774 Tau - HD 26923 

HD 26923 (V774 Tau, HIP 19859) is a “moderately active” 
GOV star that forms a loose binary system with the G8V star 
HD 26913 at an angular separation of 1' on the sky. The pro- 
j ected orb ital separation between the components is 2800 AU 
(lAbtlll988l) . Kinematically it is identified as an UMa member. 

HD 26923 is another low amplitude star. We detected peri¬ 
odicity in 15 of the 63 independent datasets with a mean V-band 
light curve amplitude of 07013. There is no previous measure¬ 
ment available for the rotation period of the star and our period 
value Pw = 11.1 d is the first one published. Some kind of mean¬ 
dering phase grouping may be perceived in the minimum epochs 
with a folding period of P = 10.7 d. However, the Kuiper statis- 
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Fig. 5. CPS results of M, A, P and 0,nin for HD 1405, HD 10008, HD 26923, and HD 29697. For M, A, and P the results from the independent 
datasets are shown with the black squares with error bars while results from the rest of the datasets are shown with points. Datasets with a constant 
brightness model (K = 0) are shown with A set to 0. The light curve minimum phases phased with Pai are shown in the fourth panel with 
black squares denoting the primary minima and grey triangles the secondary minima. This panel is shaded grey if the candidate P^] does not have 
a significant Kuiper statistic. The mean photometric period P„ is shown in the third panel as a solid line and the active longitude period P.^\ with 
the dashed line. We note that the error bars are often smaller than the plot symbols. 


tic suggests that this period is highly insignificant and does not 
allow the features to be reliably interpreted as active longitudes. 

We found a “poor” 7.0 yr cycle and looking at the mod¬ 
elled M and A results in Fig.|5] the activity variations are clearly 
quite erratic. Between 2004 and 2009 there seems to have been 
another quasiperiodic structure present in the M results with a 
period of about 2 yr. This structure is not, however, present in 
the rest of the data and was not detected by the HB method. 
iBaliunas et al.l (Il995h reported that the Ca ii H&K emission level 
of the star is variable but did not find any cycle period from their 
observations. 

5.4. V834 Tau - HD 29697 

HD 29697 (V834 Tau, HIP 2 1818) is a “very act ive” K4V star 
identified as an UMa member. [Henry et ^ d 1 995bl) found that its 
lithium abundance is consistent with the age of the Pleiades. 


We found signs of a “poor” 7.3 yr cycle in the A results of the 
star. The M results show a strong brightening trend reaching its 
tip around 2010 and turning into a slow decline thereafter. This 
might be a sign of a long cycle with a length of approximately 
30 yr or more, but claiming such a cycle at this point is highly 
premature. What is striking is the large O'T'23 amplitude seen in 
the seasonal M values. This is more than double the mean peak 
to peak value 0?0 96 of the light curve amplit ude due to rota¬ 
tional modulation. iLonez-Santiago et al.l (1201 Oh reported a pro¬ 
jected rotation velocity of vsin/ = 10.2 km s ', which implies 
a high inclination of the rotation axis. Thus the large seasonal 
variation in the mean magnitude is best explained by axisym- 
metric spot structures, e.g. a large concentration of spots around 
the pole during the brightness minimum. 

There is one apparent active longitude on the star with P^i = 
3.943 d (compare with - 3.965 d) that was particularly well 
defined before 1998. During the following decade the phase co- 


Article number, page 9 ofl24l 

















































A&A proofs: manuscript no. solan 


herence of the light curve was much poorer but after 2010 the ac¬ 
tive area has again been tightly confined around nearly the same 
rotational phase as before 1998. 

5.5. V1386 Oh-HD 41593 

HD 41593 (VI386 Ori, HIP 28954) is an “active” G9V star iden¬ 
tified as an UMa member. Our results reveal a short 3.3 yr cycle 
from its M - Ajl results. There is clear variability seen directly 
in both the M and A results but the HB method was unable to 
find any significant cycle in either one of them. The active lon¬ 
gitudes of the star with the period Pal = 8.042 d (compare with 
Pw = 8.14 d) are more easily discerned and were especially clear 
before 2006. Between 2006 and 2010 the phase coherence of the 
observed light curve minima was lost, but recently the main ac¬ 
tive longitude has reappeared at the same rotational phase where 
it was located before 2006. 

5 . 6 . V352 CMa - HD 43162 

HD 43162 (V352 CMa, HIP 29568) is the “active” G6.5V pri¬ 
mary component of a loose triple system. The other two compo¬ 
nents in the system are both M dwarfs whose proj ected distances 
from the primary are 410 AU (IChini et al.ll2014l) and 2740 AU 
dRaghavan et al.ll2010h . Age estimates of the star based on coro¬ 
nal and chromospheric activity levels indicate that it is between 
280 Myr and 460 Myr old. However, kinematically it is identi¬ 
fied as an IC member, implying that it has an age of only a few 
tens of million years. 

We found a “poor” 8.1 yr cycle in the A result s of HD 43162. 
This is reminiscent to the 11.7 yr cycle found bv iKaiatkari et al.l 
(12015h in their M results based on much of the same data as 
used in the current study. However, we could not find a pe¬ 
riod from our series of M results with FAP < 0.1 although 
there a re quasiperiodic variations. Compared to iKaiatkari et al.l 
(12015h . we used two more years of photometry in our study. The 
inclusion of extra data might have thus introduced enough insta¬ 
bility in the variations to increase the FAP above the rejection 
limit. 

We found active longitudes on the star, as did IKaiatkari et al.l 
( 2015h: the most sign ificant period P^i = 7.158 d found by 
IKaiatkari et al.l (l2015l) is much closer to the photometric ro¬ 
tation period P^ = 7.17 d than our active longitude period 
Pal = 7.132 d, and also has a higher Qk value. The active longi¬ 
tude structure connected to our Pai has been quite stable through 
the observation record with only minor back and forth phase mi¬ 
gration. 

IKaiatkari et alJ (l2Q15h also estimated the light curve period 
variation of the star using the same methods that we used. Their 
3cr range of fluctuation was Z = 0.19, which is similar to our 
value Z = 0.22. 

5 . 7 . V377 Gem - HD 63433 

HD 63433 (V377 Gem, HIP 38228) is an “active” G5V star iden¬ 
tified as an UMa member. We found two different cycles for it, a 
short 2.7 yr cycle in the M and M+Af2 and a “long” and “poor” 
8.0 yr cycle in the A and M - Ajl data. The active longitudes 
with the period Pai = 6.4641 d (compare with P^ = 6.46 d) are 
marginally detected. These are mostly visible in the data between 
2002 and 2005 and possibly also after 2012. During the first two 
observing seasons there was also a well-defined but short-lived 
drift structure seen in the light curve minima. This appears in the 


Kuiper test as a short periodicity of 6.4114 d but with (2 k = 0.16, 

which is not significant. If this drift pattern is also connected to 
active longitudes, their period variations have been quite large 
over time and the change between the two measured periods has 
been quite abrupt. 

5.8. V478 Hya - HD 70573 

HD 70573 (V478 H ya) is an “active” G6V st ar identified as a 
HLA or LA member. llsaacson & Fischeil (l2010l) estimated a very 
young age of 50 Myr for it on the basis of chromospheric emis¬ 
sion. The young age is supported by the fact that in our colour 
magnitude diagram (Fig.[T]i the star is located quite far above the 
ZAMS. There is a 6.1 Jupiter mass planeta ry companion orbitin g 
the star at a 1.76 AU semimajor axis orbit (ISetiawan et al.l200^ . 

There is a discrepancy between the values of our chromo¬ 
spheric emission index at logP^j^ = -4.488 and the p revi- 
ously reported valu es at logP^j., = -4.10 (I White et alj200’^ and 
logP^jj^ = -4.187 dlsaacson & Fischeil201ol) . which would clas- 
sify the star as “very active”. Our spectrum from February 2012 
seems to have caught the star at a state of lower-than-average 
activity. However, at the same time the star was at a brightness 
minimum implying high spottedness. It is possible that our spec¬ 
trum was simply timed so that all the major plage areas were 
located behind the limb of the star and invisible to Earth. 

There is a well-defined 6.9 yr activity cycle apparent in the 
photometry of HD 70573 that can be quite clearly traced in the 
M results. We did not And cyclicity directly in the A results, but 
variability with a slightly shorter time scale can also be seen 
in them. There has been one strong active longitude visible in 
the data throughout the observation record. This active longi¬ 
tude has been constantly migrating with respect to the average 
Pal = 3.2982 d. From the migration rates it is possible to esti¬ 
mate the seasonal coherence periods roughly as Pmigr = 3.2966 d 
between 1996 and 2000 and Pmigi = 3.2997 d both before and 
after this. Between 2010 and 2011 there was also a sudden jump 
in the active longitude phase some A<p - 0.25 forward. Each 
of these was shorter than the mean photometric rotation period 
Pw = 3.314 d. 

5.9. HD 72760 

HD 72760 (HIP 42074) is a “moderately active” KOV star iden¬ 
tified as a Hya member. There is no previous rotation period 
available for it, so our period value Pw = 9.6 d is the first one 
published. We found no evidence of any activity cycles or ac¬ 
tive longitudes on the star. In the M results there are two short 
depressions, but otherwise both the light curve mean and ampli¬ 
tude have stayed at stable levels. 

5.10. V401 Hya - HD 73350 

HD 73350 (V401 Hya, HIP 42333) is a “moderately active” G5V 
star identified as a Hya member. It had previously been thought 
to form a loose binary with the AO type star HD 73351 at an an¬ 
gular separation of 1' on the sky, but later resear ch identified the 
companion candidate as a more distant held star dRaghavan et al.l 
1201 Oh . 

There are two previous and discrepant estimates available for 
the ro tation period of H D 73350, Pmt = 6.14 d dGaidos et al.l 
I 2 OOOI 1 and Prot = 12.3 d dPetit et al.ll2008h . Using these periods 
and the known p rojected rotation velocity vsin/ = 4.0 km s“' 
dPetit et al.ll2008l) . the minimum predicted radius of the star can 
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Fig. 6. Same as Fig.|5]for HD 41593, HD 43162, HD 63433, and HD 70573. 


be calculated as either R sin i — 0.48 Rq for the shorter period 
or Rsini = 0.97 Rq for the longer period. Considering that 
the spectral type is close to solar makes the longer period more 
likely. More conclusive proof of the longer period value comes 
from our preliminary TSPA analysis of the photometry which re¬ 
vealed the roughly 12 d period as the primary periodicity in the 
data and the roughly 6 d period as its first overtone, which re¬ 
sults from the secondary minima often present in the light curve. 
Our final mean photometric rota tion period from the CPS was 
Pw = 12.1 d. It is remarkable that iPetit et al.l (l2008h were able to 
find this period from only 13 nights of data by using the period 
as a fitting parameter in their Zeeman Doppler imaging. 


We found a short 3.5 yr cycle in the mean magnitudes M. The 
cycle looks fairly stable and although we did not find it in the A 
results, there is similar behaviour visible in them as well. The 
light curve amplitude is typically higher during the times when 
the mean brightness is at its minima. A candidate folding period 
at P = 12.59 d can be found by the Kuiper method and corre¬ 
sponds to a phase structure seen before the year 2002. However, 
this period has a high Qk value and cannot be identified with 
stable well-defined active longitudes. 

5.11. DX Leo-HD 82443 


iPetit et aP (l2008b also tried to estimate the absolute value of 
the surface differential rotation by similar means but obtained 
a value with large error bars, AQ = 0.2 + 0.2 rad d *. If we 
interpret the range of our rotation period fluctuations, Z = 0.40, 
as a direct measure of the differential rotation coefficient k, we 
obtain AQ = 0.21 rad d ' from our measu rements. A l thoug h this 
is a high value, it fits the mean estimate of iPetit et al.l (l2008h very 
closely. 


HD 82443 (DX Leo, HIP 46843) is a thoroughly studied “ac¬ 
tive” KIV star. It forms a loose bi nary system with a M4.5 
type dwarf at a distance of 1716 AU dPoveda et al.lll994ll . Kine¬ 
matically it has been grouped as a THA, Pie, or LA member 
and its lithium abun dance is consistent with the Pleiades age 
(iMontes et alJl200T^ . 


A s hort activity cycle ha s been found by many authors for 
the star. lBaliunas et al.1 (1 1995) reported a 2.8 yr c ycle in the chro¬ 
mospheric emission while Messina et al.l (Il999h reported a cycle 


Article number, page 11 ofl24l 




























































A&A proofs: manuscript no. solan 


HD 72760 



year 


V401 Hya - HD 73350 



year 


DX Leo - HD 82443 



year 



year 


Fig. 7. Same as Fig.|5]for HD 72760, HD 73350, HD 82443, and HD 82558. 


of 3.89 yr and iMessina & GuinanI (l2002l) of 3.21 yr, both from 
photometry. Our results for the variation of M and A show at 
first glance quite erratic behaviour, but a “fair” 4.1 yr cycle is 
still firmly present in both of them in addition to the M - Af2 
and M + Af2 results. Furthermore, we were able to find a new 
“long” 20.0 yr cycle from the M values. 

The active longitude behaviour of HD 82443 is both well 
developed and complex. First, between the 1993 and 1994 ob¬ 
serving seasons there was a flip-flop event on the star. While the 
event was underway, both of the active longitudes were present 
for a while but one of them gradually weakened and eventually 
disappeared entirely, while the other one took over as the new 
primary active longitude. Later around 1996, this new active lon¬ 
gitude started to migrate with respect to the average active lon¬ 
gitude period Pal = 5.4147 d with a period of Pmigr = 5.4243 d 
estimated from the migration rate. This migration pattern even¬ 
tually stopped and turned into an opposite migration between 
2003 and 2005 with a period of Pmigr = 5.4029 d. Each of these 
periods is shorter than or equal to the mean photometric period 
Pw = 5.424 d. Since 2005 the star has shown one stable active 
longitude with some migration back and forth but mostly staying 


close to a single phase in the average Pai = 5.4147 d rotational 
frame of reference. 

IMessina et al.l (1 19991) estimated the surface differential rota¬ 
tion of the star from obs erved variations in the pho tometric pe¬ 
riod as AP/P > 0.04 and IMessina & GuinanI (l2003l) similarly as 
AP/P = 0.024. These estimations are in line with our small value 
Z - 0 .054 calculated using similar means. IMessina & GuinanI 
(l2003l) also correlated their estimated period values with the 
phase of the activity cycle and claimed that the differential ro¬ 
tation would be antisolar. No trace of such behaviour can be 
confirmed by our results and it remains impossible to specify 
the sign of the surface differential rotation. 


5.12. LQ Hya-HD 82558 


HD 82558 (LQ Hya, HIP 46816) is another “very ac¬ 
tive” KOV star with a rich reference history. It ex- 
hibits strong chromospheric emission modulated by rotation 


mbits strong cnromospneric emission moaulatea by rotation 

( Fekel et a l.fl986bl: Shassmeier et aI.II993HMontes et alJ2001at 


I Alekseev & Kozlo^ l2002l) in such a way that stronger emis¬ 

sion is concentrated on the same phases as the dark photospheric 
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spots (iFrasca et alJ 120081 : ICao & Gull20l4 iFlores Soriano et aP 
l2015h . It has also been observ ed to have strong flares in the op¬ 
tical and X-ray wavelengths (iMontes et alJ 119991: ICovino et alJ 
l2Q0lh . The star is very young as demonstrated by the identifica¬ 
tion as an IC membe r. The high lithium abundance also points to 
a Pleiades-type age (iFekel et al.lll986bl) . 

Comple x surface magnetic fi elds were directly observed on 
the star by iDonati et alJ (Il997h . Based on Doppler imaging 
and Zeeman Doppler imaging, the spot activity appears to be 
concentrated in two areas, one at low latitudes and the other 
near the pole ( Strassmeier et al.lll993l: Rice & Strass meieill 19981: 


Donatilll999HDonati et al.ll2003bl: iKovari et ahl 2004 ilCole et alJ 


2015lf Comparison of spot longitudes bv ICole et al. (l20L5h 
tween Doppler imagi ng re sults and the p hotom etric results by 
iLehtinen et al.l (1201 2h and lOlspert et al.l (1201 5h revealed good 
agreement between the two observational methods. 

Activity cycles have been widely reported for HD 82558 
with lengths mostly grouping into three distinct ranges. Cycle 
lengths repo rted previously from th e photometry are: 6.24 yr by 


Jetsul (Il993l): 6.8 vr and 1 1.4 yr bv lOlah et al.l (20 001): 3.4 yr by 


Olah & Strassmeieij(l2002h : 7.7 yr an d 15 yrby Berdvugina et al.l 
( 2002h : 3.2 vr. 6.2 vr. and 11.4 yr W Messina & Guinan (l2002l) : 


3.7 yr, 6.9 yr, and 13.8 yr by iKovari et al.l (12004 : 2.5 


yF 

3.6 yr, and a c ycle in creasing in lengt h from 7 yr to 12.4 y r 
bv lOlah et al.l (l2()09l): and 13 yr by ILehtinen et al.l (l2012h . 
iBerdvugina et al.l (l2002h also reported a 5.2 yr cycle from the 
flip-flop behaviour of the active longitudes, though this does not 
seem likely in the light of our results. 

We were able to find a highly significant “long” cycle in our 
data but its exact length varied quite a bit depending on which 
set of photometric results we looked at. The cycle lengths re¬ 
trieved by the HB method were 17.4 yr for the M results, 14.5 yr 
for the A results, 15.8 yr for the M - Ajl results and 18.0 for 
the M + A/2 results. These doubtless represent the same under¬ 
lying cycle but its long duration and seemingly non-stationary 
behaviour mean that there is considerable instability in the time 
series analysis. Our results also point to a longer cyc le than any 
of the past results. While the analysis of Lletsul (119931) revealed a 
good 6.24 yr fit to the light curve mean magnitude between 1984 
and 1992, the long-term variations of our M and A results have 
shown increasingly sluggish behaviour after that. Since 2000 the 
average mean brightness of the star has been continuously in- 
creasing with no sign of turning back. In this light the result of 
lOlah et akl (|2009|) of a cycle with an increasing length seems like 
an apt description of what is happening on the star. We also note 
that our M results show small amplitude oscillations with a time 
scale of 2 to 3 years similar to of the shortest cycles reported by 
other authors. The HB method was, however, unable to detect 
these shorter cycles from our data. 

HD 82558 shows only occasional active longitude behaviour. 
Simil arly to our previous s tudy using only the V-band photom¬ 
etry (ILehtinen et al.l 1201^ . we were able to find coherent ac¬ 
tive longitudes between 2003 and 2008 but for the rest of the 
time no other active longitude patterns could be found. A short¬ 
lived pattern seen for some years before 1995 at a period of 
P - 1.68929 d or F = 1.61208 d could be found from our previ¬ 
ous data but is absent from our current results. The inclusion 
of more data also meant that the significance of the detected 
period at Fai = 1.603733 d (compare with Fw = 1.6044 d) 
dropped quite a bit from our previous study. The la c k of s table 
active longitudes contra sts with the results ofLletsul (Il993l) and 
IBerdvugina et al. d200^ but agrees with iDonati et al.l (1200^ 
and lOlspert et al.l (1201 5l) who find n o persistent lo n gitudi nal ac¬ 
tivity concentrations. Furthermore, lOlsnert et al.l (l20L5h . who 


analysed the same V-band data we did, found that no consis¬ 
tent light curve period could be identified from the data using 
correlation times longer than approximately 230 d. 

Differential rotation of HD 82558 has been estimated from 
the de tected photometric period variations as Z = 0.015 (iJetsul 
Il993h and Z = 0.020 (ILehtinen et al.l 1201^ using the same 
methodology we usre in this paper and as AF/F_= 0.013 
(iMessina & Guinanll2003l) and AF/F = 0.025 (lYoull200^ from 
the total observed period range. These are all in line with our 
current estimate of Z = 0.017. Other methods have produced 
smaller differential rotation estimates. [Berdvugina et al.l (l2002h 
found the value k - 0.002 from tra cing the migr a tion r ates of 
their detected active longitudes and IKovari et al.l (l2004l) got a 
value of k - 0.0056 b y cross-corre l ating th eir Doppler imag¬ 
ing temperature maps. IDonati et al.l (12003 al) used surface dif¬ 
ferential rotation as an optimized parameter in their Zeeman 
Doppler imaging and found values of 0.0037 < k < 0.049 for 
the Stokes I inversion and -0.013 < k < 0.051 for the Stokes 
V inversion. They concluded that this range of values possi- 
bly represents physical ch anges in the differential rotation rate. 
IMessina & GuinanI (l2003l) correlated their local period estimates 
with the phase of their 6.2 yr cycle and claimed antisolar differ¬ 
ential rotation. As in the case of HD 82443, this behaviour is not 
seen in our results. In any case, it appears safe to say that the 
differential rotation of HD 82558 is small and the star exhibits 
close to rigid body rotation. 

5.13. NQ UMa - HD 116956 

HD 116956 (NQ UMa, HIP 655 15) is an “active” G9y star iden¬ 
tified as a TWA or LA member. ILehtinen et al.l (1201 ll) suggested 
the existence of a 3.3 yr cycle in a subset of the same V-band 
photometry as in the present study. We confirmed this here af¬ 
ter finding a cycle of 2.9 yr in our M and M - A/2 results. We 
also found an additional “long” cycle of 14.7 yr in the A results. 
The same two active longitudes with the period Fai = 7.842 0 d 
(compare with F* = 7.86 d) found bv ILehtinen et al.l (1201 ll) are 
also clearly present in our results and span the whole observation 
record without interruption. The active longitudes have under¬ 
gone two phase jumps smaller than A0 = 0.25 and their order of 
strength has switched around a few times, but for the most part 
the active longitudes are very constant. We previ ously estimated 
the d ifferential rotation of the star as Z = 0.11 (ILehtinen et al] 
1201 ll) which is identical to the value we found from our new 
data. 


5.14. KULib-HD 128987 

HD 128987 (KU Lib, HIP 71743) is a “moderately active” G8V 
star with a low light curve amplitude. Kinematically it is iden¬ 
tified as an IC member. We found a 5.4 yr cycle from its M, 
M - A/2, and M + A12 results but did not find evidence of active 
longitudes. 

5.15. HP Boo-HD 130948 

HD 130948 (HP Boo, HIP 72567) is a “moderately active” 
F9IV-V star. It is orbited by a binary s ystem of two L-t ype 
brown dwarfs at a projected distance of dPotter et'ani2002l) . It 
has not be en identified as a memb er of any 48 AU kinematic 
group and iMaldonado et al.l (1201 Ol) simply labelled it a young 
disk object. Estimates based on chromospheric and coronal ac¬ 
tivity imply ages between 190 Myr and 870 Myr. The spec- 
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Fig. 8. Same as Fig.|5]for HD 116956, HD 128987, 130948, and HD 135599. 


tral classification of the star would allow it to be classified as 
a more evolved subgiant. We calculated its minimum predicted 
radius from the rotation period and the projected rotational ve¬ 
locity vsini = 8.54 km s“' (iMartmez-Arnaiz et al.l 1201^ as 
R sin i - 1 3Rq. Based on this value, the age estimates, and be¬ 
cause it falls directly on the ZAMS in our colour magnitude dia¬ 
gram (Fig.[TJ, it is most probable that this star is a main-sequence 
object. 

We found a “poor” 3.9 yr cycle from the M and M + AI2 
results of the star. There is no evidence of any presence of active 
longitudes. 


5.16. V379 Ser - HD 135599 


HD 135599 (V379 Ser, HIP 74702) is an “active” KOV st ar. It 
was identified as an UMa member by Gaidos et al.l (l2000t) but 
this was refuted bv iMaldonado et aTI (1201 Ol) who described it as a 
“probable non-member”. Age estimates based on chromospheric 
and coronal activity range between 200 Myr and 1340 Myr sug¬ 
gesting an older age for the star. We found a long 14.6 yr cycle 
from our M and M + Ajl results. There is a suggestive correla¬ 


tion between the M and A results since the light curve amplitude 
reached its maximum values around the same time as the mean 
brightness dipped to its minimum. Still, we did not find any cy¬ 
cle from the A results, nor did we find any evidence of active 
longitudes. 

5.17. V382 Ser-HD 141272 

HD 141272 (V382 Ser, HIP 77408) is the “moderately active” 
G9V primary component of a loose binary star. The secondary 
component of the system is an M dwarf at a proje cted distance 
of 350 AU from the primary (lEisenbeiss et alj2007h . The system 
is identified as a member of HLA or LA. We were able to find a 
“poor” 6.4 yr cycle from our M results but did not find evidence 
of active longitudes. 

5.18. V889 Her-HD 171488 

HD 171488 (V889 Her, HIP 91043) is a “very active” G2V star 
with a ri ch reference histo r y. It is among the youngest stars in our 
sample. IStrassmeier et al.l (l2003h estimated an age of 30 Myr by 
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Fig. 9. Same as Fig.|5]for HD 141272, HD 171488, HD 180161, and HD 220182. 


comparing it with isochrones and 50 Myr from the lithium abun¬ 
dance, while iFrasca et al.l (l2010h estimated an age of 50 Myr 
from isochrones. Kinematically the star is identified as a LA 
member. 

The chromospheric emission of HD 171488 has been ob¬ 
served to_be_concenttatedon_the_same phases as the dark 
spots (iMulliss & Bonnl [1994 iFrasca et^ l2010ll . Doppler im¬ 


ages have been published by a number of authors and they 
all indicate that the spot activity is dominated by one or a 


few Dolar snots (IStrassmeier et al. 20031 Marsden et al. 

2006 

Jarvinen et al. 2008; iJeffers & Donatil 2008; 

Huber et al 

2009 

Frasca et al.||2010|l. 


There is a clearly defined 9.5 yr cycle present in the M, 
M - Af2 and M + A/2 results of the s tar. This had previously 
been identified hv Llarvinen et al.l (l2008l) as a 9 yr cycle. Two full 
cycles of roughly the same length can be seen in the M results 
superimposed on a slight dimming trend. Variations can also be 
seen in the A results, including two fast changes in 1995 and 
2014, but no cycle could be found from them. There is also a 
pattern seen in the long-term evolution of the P results but we 
are not able to claim anything about its nature. 


An active longitude with the period Pal = 1.33692 d (com¬ 
pare with Pw = 1.345 d) has been present on the star between 
1994 and 1998 and later between 2003 and 2008 with minor 
migration back and forth. For the rest of the time the phase dis¬ 
tribution of the light curve minima has been less coherent. Ten¬ 
tatively, the active longitude appears to have been most clearly 
visible during the periods of increasing mean brightness or de¬ 
creasing spottedness, but it is still too early to claim this with 
certainty. 

The absolute value of surface shear on HD 171488 due to 
differential rotation has been estimated from Zeeman Doppler 
imaging by usin g the differeri t ml ro tation rate as an opti¬ 
mized parameter. Marsden et al.l (l2006h found a value of AQ = 
0.40 r ad d ' from their Stokes 1 inversion and Lleffers & Don^ 
(l2008h the values AQ = 0.52 rad d * from their Stokes I in¬ 
version and AQ = 0.47 rad d ’ from their Stokes V inversion, 
which correspond to differential rotation coefficients k - 0.084, 
k = 0.11, and k - 0.10, respectively. They are somewhat larger 
than what our value of Z = 0.054 suggests, possibly implying a 
narrow latitude range for the spots on the star. 
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Fig. 10. Same as Fig.|5]for SAO 51891. 


5.19. MVDra-HD 180161 


HD 180161 (MV Dra, HIP 94346) is a “moderately active” G8V 
star identified as a H ya member. It s rotati on period was estimated 
as Prnt = 9.7 d b y ^^dos et al.l (l2000l) . but Prot = 5.49 d by 
IStrassmeier et al.l (120001). Using the protected rotational velocity 
vsin/ = 2.20 km s~* dMishenina et al.l ^12h . we can calculate 
the minimum expected radius of the star as R sin i — Q.42 Rq from 
the longer period and as R sin i = 0.24^0 from the shorter period. 
This strongly favours the longer period, which was confirmed by 
our TSPA analysis. The final mean photometric rotation period 
from the CPS analysis was = 9.9 d. 


Some erratic short-term variation can be seen in the M results 
of HD 180161 but no cycle could be found from the star. Neither 
did we find any evidence of active longitudes. 


5.20. V453 And-HD 220182 

HD 220182 (V453 And, HIP 115331) is an “active” G9V star. It 
has not been assigned to any kinematic group, but gyrochronol- 
ogy and the chromospheric and coronal activity point to an age 
of more than 100 Myr. The full range of published age estimates 
is from 40 Myr to 550 Myr. 

We found a “long” cycle of 13.7 yr from the M and M H-A/2 
results. By visual inspection a roughly 6 yr cycle might also be 
present in the M results but we did not find any trace of it with 
the HB method. A persistent active longitude has been present 
on the star for the full span of the observing record, possibly 
excluding the last few observing seasons where the data was 
sparser. There have been two major migration patterns around 
the average active longitude period P^\ - 7.6200 d. From the 
migration rates we estimated that between 2000 and 2002 the 
active longitude followed a shorter rotation period of approxi¬ 
mately Ptnigr = 7.610 d and both before and after that a longer 
period of approximately Pmigr = 7.625 d, which have all been 
shorter than the mean photometric period Pw = 7.68 d. A possi¬ 
ble secondary active longitude separated by A0 = 0.5 from the 
primary might also be present, but it is not particularly strong. 


5.21. V383 Lac-SA0 51891 


SAO 51891 (V383 Lac) is an “active” KIV star, although ac¬ 
cording to the previous ly estimated chrom ospheric emission in¬ 
dex logPjjj. = -4.05 (IWhite et alJl2007 h i t coul d also be clas¬ 


sified as a “very active” star. Biazzo et al.l (l2009h observed that 


the Ca II H&K, Ca ii IRT and He emission is concentrated on the 
same rotational phases as the photospheric spots. They also ob¬ 
served significant variation in the Ha emission, but did not see 
any dependence between it and the rotational phase. They rea¬ 
soned that the Ha emission might be dominated by processes 
like microflaring causing intrinsic variation and masking away 
the rotational modulation. T he star is identified as a LA member 
and iMulliss & Bot)^(ll994l) concluded that its lithium abundance 
is consistent with an age younger than the Pleiades. 

Our photometric time series of SAO 51891 is dominated by 
the long gap separating the observations gathered between 1994 
and 1996 from the observations gathered from 2011 onwards. 
Consequently, the CPS results consist of only 15 independent 
datasets with reliable model fits and we were not able to find ac¬ 
tivity cycles or active longitudes from them. Nevertheless, there 
is evidence of a possibly quite long cycle seen in the M results. 
The mean V-band magnitude observed around 1995 was approx¬ 
imately 0.2 magnitudes dimmer than that observed during the re¬ 
cent years and its trends were opposite during the two parts of the 
photometry. The amplitude of the M variations is exceptionally 
large compared to the rest of our sample stars. It is thus expected 
that with continued observations a distinct large amplitude cycle 
will become evident on SAO 51891. Likewise, the best period 
P - 2.40 d found by the Kuiper method reveals a tightly con¬ 
fined feature from the light curve minima during the first three 
years of observations. It is thus also possible that well-defined 
active longitudes will surface from the star as more observations 
are gathered. 


6. Discussion 

6.1. Differential rotation 

In Sect. 14.11 we touched on the problems of interpreting the Z 
values as estimates of differential rotation. To quantify the ef¬ 
fect of spurious fluctuations in the rotation period estimates, we 
simulated a set of periodic test time series with constant periods 
and using a range of values for their length, the number of data 
points contained in them, and the level of noise in the data. We 
then performed our period search for these simulated time se¬ 
ries and calculated the resulting values of Z. It is evident that the 
three varied parameters (the number of data points per dataset 
«data, the number of complete rotations within each dataset 11 ^ 1 , 
and the ratio of noise to the signal amplitude e) contribute to the 
instability of the period determination. We chose their test ranges 
to be similar to what is expected for our ground-based automated 
photometry: 10 < iidata ^ 100, 2 < iiiot < 20, and 0.02 < e < 0.5. 

We found that the spurious fluctuations Zjpu quite expectedly 
increased as the data points became scarcer, the noise level in¬ 
creased, or the datasets included fewer complete rotations. There 
is a linear relation of Zspu to each value of and e. We 

found that the best fit to the simulation results was obtained by 

Zspu = 16.0 en„'t(«da‘ta + 0-023) (10) 

The fit is demonstrated in Fig. [TT] for a selection of computed 
Zspu values. 

If we calculate the mean values of iidata and «rot per dataset 
and the mean ratio e of observational errors to the light curve 
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Fig. 11. Example slices of the fit of Eq. [T0]onto a selection of Zspu values. Left: Zsp^ computed for Udata = 10 and e € {0.05,0.2,0.5) ranging from 
the darkest to the lightest grey. Centre: Z^pi, for = 2 and e e {0.05,0.2,0.5) from the darkest to the lightest grey. Right: Zjp^, for = 2 and 
Mdata 6 {50,20,10) from the darkest to the lightest grey. The scale in Zspu is identical for all the panels. 


amplitude for the CPS results of our stars, we find that the pre¬ 
dicted values of Zspu are always in the range of a few percent of 
the Z values calculated from observations. The Z estimates seem 
thus fairly rigid against numerical instability and correcting them 
for such small effects is not reasonable. There is also reason to 
believe that within the time scale of a few rotations, as used 
for our datasets, active region growth and decay also produce a 
relativ ely minor effect on the period fluctuations dPobson et alJ 
Il990l) . We can thus say that the uncertainty in relating the Z from 
our analysis to the relative differential rotation coefficient k is 
dominated by the unknown latitude extent of the spot areas on 
the stars. The exact values of k may remain unknown, but we 
can use Z as a proxy measurement by using the proportionality 
Z cc k. 

Despite the uncertainties associated with measuring the dif¬ 
ferential rotation, there is a clear dependence between Z and 
the rotation period. Stars that have shorter rotation periods have 
smaller values of Z and thus appear to have smaller differential 
rotation. For our stars we find a linear relation in the logarithmic 
scale as 

logZ =-1.85 H-1.36 logPpot- (11) 

The stars and the fit are shown in Fig. [12] as the squares and the 
solid line. 

This relation can be compared to previous results obtained 
by other authors. iHenrv et al.l (Il995ah performed a study where 
they traced the rotation of light curve features identified as dis¬ 
tinct spots. They estimated k using the difference AP between 
the largest and smalle st detected p eriods. For 87 stars taken from 
their study and that of fH^ (1199 Ih they found the relation 

logk = -2.12 H- 0.76 logPpot - 0.57 F, (12) 

where F is the Roche lobe filling factor. iDonahue et al.l (Il996h 
investigated rotational variation seen in the Ca ii H&K line emis¬ 
sion on 36 active stars and the Sun and estimated differential 
rotation by finding AP from seasonally computed periodograms. 
They found a scaling law AP oc P^^^, corresponding to scaling 
the differential rotation coefficient as k oc AP/P oc 

The relat i ons f ound by iHenrv et^ (Il995ah and 
iDonahue et alJ d 19961) are both shown in Fig. [T2| in com¬ 
parison to our results. The dashed lin e and grey points denote 
the relation by IDonahue et al.l (Il996|) and the s t ars the y used. 
The dotted line marks the relation by lHenry et alJ d 1995 ah within 


the shown period range. Since all our stars lack close binary 
companions, we haye set F = 0 for Eq.[T2| The relation predicts 
smaller yalues of diff erential rotation com pared to the other re¬ 
sults perhaps because iHenry et~^ d 1995 ah traced the migration 
patterns of indiyidual longer-liyed spots instead of looking for 
the full period yariation, thus potentially downplaying the effect 
of short-liyed spots. 

A further study was made by iBarnes et alJ d2005h based on 
Doppler imaging who found the scaling of the absolute equator 
to pole shear AD to the rotation rate Q to be AQ cc D° This 
corresponds to k oc which is again qualitatiyely similar to 
results from the other studies in that it shows increasing differ¬ 
ential rotation towards slower rotators. In general we find that if 
we haye k oc P^^^ for the differential rotation coefficient, we get 
AQ oc Q’' = Q^“^ for the absolute shear. The power law indices 
H and V are listed for all the mentioned studies in Table |7] 

The yalues of fi and v show considerable scatter between the 
different studies. The differences between the exact obtained yal¬ 
ues may originate from different choices of data analysis strate¬ 
gies and the physical rotation proxies being tracked, and from 
simple random scatte r due to limited sample sizes. In the case of 
our sample and that of iDonahue et alJ (Il996l) . eyen the power law 
nature of the underlying scaling cannot be fully yerified since 
they both coyer only approximately one dex in Frot- 

Neyertheless, at least in the case of p, the results can be 
interpreted to qualitatiyely agree with each other since they 
show slower rotating stars to haye stronger differential rota¬ 
tion. In the case of the dependence of AQ on Q the results 
are less clear, but as they are scattered closer to v = 0 they 
may still suggest an agreement w ith theoretical results such as 
those of iKiiker & Rudigeil (1201 ll) who found AQ to haye lit¬ 
tle dependence on the rotation rate in t heir mean field models 
of l ower main-sequenc e stars. Recently iReinhold et al.l (1201 3l) 
and [Reinhold & Gizoiil (12015l) performed a massiye study of 
24 124 spotted Kepler field stars and measured differential ro¬ 
tation from two simultaneously obseryed rotation periods. Their 
results agree quite well with the theoretically predicted indepen¬ 
dence of AQ from Q and ar e yery close to the obseryational re¬ 
sults of lHenry et al.l (Il995al) . 

Ac cording to iBarnes et al.l (l2005l) and ICollier CameronI 
(l2007l) there is a strong dependence of AQ on the temperature 
of the star. We were unable to find eyidence for this in our re¬ 
sults, only finding a Pearson correlation coefficient r - -0.21 
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Fig. 12. log Z vs. log Plot for our sample stars shown with the squares. 
The colour coding denotes the B — V colour of the stars. The fit to our 
results ('Ea .Il lb is shown with th e solid line. For comparison we plot the 
stars from iDonahue et ^ ( Il996h as grey points and the fit t o these stars 
with th e dashed line. We also show the fit to the stars of iHenrv et al.l 
( Il995ah wi th the dotted line for the log P^ot range covered by our work 
and that of iDonahue et alJ lll996h . 

Table 7. Power law indices p and v for differential rotation fits k oc 
and An oc £2'' from the current study and the literature. 


p 

V 

Reference 

1.36 

-0.36 

this work 

0.76 

0.24 

Henrv et al. (1995a) 

0.3 

0.7 

Donahue et al. (1996) 

0.85 

0.15 

Barnes et al. (2005) 

0.71 

0.29 

Reinhold & Gizon (201^ 


between the B — V colour and Z/Prot AQ. This lack of cor¬ 
relation can also be seen from the colour co ding for B - V in 
Fig. [121 Our results agree better with those of iKuker & Riidig^ 
(|2Q1 ill , who found only weak temperature dependence of AQ 
for stars with temperatures within the range 3500-6000 K and 
a steeper dependence only for stars hotter than this. The two- 
staged te mperature dependence was also observationally con¬ 
firmed by iReinhold et al.l d2013ll whose sample included stars 
from both sides of the 6000 K divide. 

6.2. Active longitudes 

We were able to identify reliable active longitudes for roughly 
half of our stars. The values found for P^i are shown in Fig. [T3] 
in comparison to the weighted mean photometric periods P^. 
All the periods are normalized in the plot by dividing them by 
Pw The heights of the boxes corresponding to the period values 
denote the +lcr error limits of the period estimates. 

The Pal values are smaller than the P^ values nearly without 
exception. In the case of HD 41593 the relative difference be¬ 
tween the two periods is as large as 1.1% which corresponds to 
a lap time of 1.9 yr between Pai and P^. This raises the question 
of whether there is an underlying difference between the coher¬ 
ence period of the active longitudes and the photometric rotation 
period of the stars. This situation is possible if the active longi¬ 
tudes represent an activity generating structure rotating with a 
different angular velocity compared to the stellar photosphere. If 


Fig. 13. Comparison of the active longitude periods Pai (dark boxes) to 
the photometric mean rotation periods P„ (light boxes), both normal¬ 
ized to the values of P*. The heights of the boxes show the estimated 
±1(T error bars of the periods. 


individual spots or spot groups are formed at the active longi¬ 
tudes but start to follow the photospheric rotation once formed, 
we would expect there to be separate signals in the photometry 
for the periods Pai and P^. 

To get an idea whether there is any significant difference be¬ 
tween Pal and Pw, we modelled the error distributions of P* as 
Gaussians )V(Pw,crp^) with their standard deviations as the P^ 
standard errors of the mean. We then calculated the two sided 
tail-area probabilities p - Pr(|P-Pw| > |Pai-Pwl) that a random 
sample P drawn from such a distribution has a larger deviation 
from Pw than Pai has. This gives us a measure analogous to the 
Qk and FAP used in the Kuiper test and the HB method. The 
process ignores the fact that the Pai values also have error distri¬ 
butions associated with them, but since they are always narrower 
than CTp the simple p values still provide a useful heuristic to 
assess the distinctiveness of Pai from P„. 

The relative differences of Pai from P*, 6P - Pai/Pw - 1, 
and the tail area probabilities p for each found active longitude 
period are listed in Table|8] In the case of HD 43162, HD 63433, 
HD 82558, and HD 116956 there is little evidence of observed 
difference between Pai and P^. On the other hand, in the case of 
HD 1405, HD 29697, HD 41593, HD 70573, and HD 171488 the 
directly detected values of Pai have p < 0.01 and for two more 
stars (HD 82443 and HD 220182) p < 0.05. For these stars, a 
single underlying period appears inadequate to describe both Pai 
and Pw at the same time. 

The observed differences between Pai and P„ can have sev¬ 
eral different interpretations. It might be that the magnetic struc¬ 
tures that sustain the active longitudes reside at deeper levels in 
the stellar interior than the formed starspots and trace the rotation 
at those depths. In this case the difference between Pai and P* 
would reflect radial shear in the outer layers of the stellar inte¬ 
rior. Helioseismological results have shown that in the Sun there 
is a negative gradient of angular veloci ty on all latitudes just be¬ 
low the photosphere (iHowe et al.ll2000l) . This behaviour of inter¬ 
nal rotation would predict shorter rotation periods for structures 
rooted at deeper levels, just as observed for the active longitudes 
on our stars. 

On the other hand, there is theoretical evidence from direct 
numerical simulations that a stellar dynamo may operate with 
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a non-axisymmetric mode that propagates longitud inally with 
respe ct to the rotational reference frame of the star (ICole et alJ 
|2014|) . If this is the case for our sample stars, what we are seeing 
as the migrating active longitudes are signs of azimuthal dynamo 
waves. Such a dynamo wave would feed a certain longitude in 
the rotational frame of P^i with spot structures, which could then 
start to drift with the local surface rotation. These are the az¬ 
imuthal counterparts to the latitudinal dynamo waves seen on 
the Sun and ar e responsible for the well-known butterfly diagram 
(iRadledI 19861) . 

Similar drift patterns have b e en ob served before from ac¬ 
tive giants stars. [Hackman et'aP (1201 Ih observed that the ac¬ 
tive longitude on II Peg has rotated with a shorter period than 
the stellar surface, as is the case for our stars, while FK Com 
(iHackman et al.ll^l3h seems to have followed an opposite pat¬ 
tern of the active longitude rotating more slowly than the stellar 
surface. What still need to be properly explained are, on the one 
hand, the observed values of the relative period differences 6P 
between Pal and P^ and, on the other, the variable migration 
trends of the active longitudes seen on some stars. In Fig. [13] 
we have sorted the stars in the order of increasing P^ and no 
immediate relation between rotation period and the 6P value is 
apparent. 

We note that it does not appear justified to interpret the differ¬ 
ences seen between the periods Pai and Pw as measures of d iffer- 
ential rotation, similarly to how iReinhold & Gizoril (l20L5l) used 
two simultaneous photometric periods detected in single quad¬ 
rants of Kepler data for their differential rotation estimation. If 
such an interpretation were true for our much longer time series, 
often showing very stable active longitudes, we should expect a 
corresponding bimodality to be visible in the distribution of the 
local P estimates produced by the CPS analysis. A look at the 
estimated P values in Figs. BI-fTOl reveals instead only unimodal 
distributions centred at the Pw values. 

A second significant result concerning the active longitudes 
is their presence on only the more active stars. This is shown in 
the right-hand panel in Fig.[T4| where we have circled those stars 
with activity cycles that also have reliably detected active longi¬ 
tudes. There is a sharp divide at approximately logPj^j, = -4.46 
between active stars showing active longitudes and less active 
stars with poor or no evidence for them. In our sample there are 
only two exceptions to this pattern. One of them is HD 70573, 
which has a chromospheric emission level slightly below the 
above-mentioned divide and for which we have detected active 
longitudes. The other exception is SAO 51891, which shows 
high chromospheric activity at logPjjj. = -4.327 but for which 
we could not determine any active longitudes or definite activity 
cycles. We note that the lack of these detections is likely due to 
the long time gap in the observations from the star. 

These results suggest that there are two well-defined do¬ 
mains of magnetic held geometry present on the late-type ac¬ 
tive stars that are defined by their activity level. On the less ac¬ 
tive stars the large-scale magnetic held appears to prefer an ax- 
isymmetric configuration, whereas on the more active stars the 
held configuration is dominated by a non-axisymmetric mode. 
This observation can be related to the numerical results of 
iTuominen et al.l (Il999h who found their mean held dynamo so¬ 
lutions to switch from axisymmetric to non-axisymmetric with 
increasing Taylor numbers, i.e. higher rotation rates. 

6.3. Activity cycle lengths 

With a sample of estimated cycle periods it is possible to 
study the connection between the cycle lengths and other stel¬ 


Table 8. Relative differences 6P between the active longitude and pho¬ 
tometric periods and the tail-area probabilities p that the underlying 
photospheric rotation periods have a larger deviation from the estimated 
Rw than Rai- Values of 6P <0 indicate P.jx < P^. 


Star 

Pal [d] 

5P 

P 

HD 1405 

1.752212 

-0.23% 

O 

o 

1 

HD 29697 

3.9433 

-0.55% 

2.2 ■ lO-'* 

HD 41593 

8.0417 

-1.1% 

1.9 - 10-3 

HD 43162 

7.1323 

-0.50% 

0.35 

HD 63433 

6.46414 

0.033% 

0.96 

HD 70573 

3.29824 

-0.49% 

2.3 ■ 10-® 

HD 82443 

5.41471 

-0.18% 

2.4 ■ 10-3 

HD 82558 

1.6037330 

-0.038% 

0.14 

HD 116956 

7.84203 

-0.23% 

0.26 

HD 171488 

1.336923 

-0.63% 

7.0- 10-" 

HD 220182 

7.62002 

-0.76% 

1.2- 10-3 


Table 9. Semi-empirical convective turnover times and inverse Rossby 
numbers. 


Star 

Tcld] 

logRo ‘ 

HD 1405 

22.53 

2.21 

HD 10008 

19.27 

1.55 

HD 26923 

7.45 

0.93 

HD 29697 

23.71 

1.88 

HD 41593 

19.85 

1.49 

HD 43162 

15.51 

1.43 

HD 63433 

13.82 

1.43 

HD 70573 

10.48 

1.60 

HD 72760 

19.05 

1.40 

HD 73350 

12.28 

1.10 

HD 82443 

18.58 

1.63 

HD 82558 

22.36 

2.24 

HD 116956 

19.51 

1.49 

HD 128987 

15.35 

1.31 

HD 130948 

7.77 

1.09 

HD 135599 

20.34 

1.66 

HD 141272 

19.41 

1.25 

HD 171488 

10.13 

1.98 

HD 180161 

19.51 

1.39 

HD 220182 

19.41 

1.50 

SAO 51891 

20.80 

2.03 


lar para meters. This was done by B^randenburg et al.l (Il998h and 
later bv ISaar & Brandenburg! (1 19991) who studied the relation of 
Frot/Fcyc to the inverse Rossby number Ro ' and logRjjj^. They 
were able to group their stars into various activity branches char¬ 
acterized by the different behaviour of Piot/Fcyc and suggested 
that these branches are related to a dynamo evol ution sequence. 


Slightly different types of studies were done byjBaliunas_e^^ 


(1996™, w ho compared Pew. IP mi to l/Pmt, and iBohm-Vitense 


(2007) and lOlah et al.l (l2009t) . who also compared Peye directly 
to Prot- In our data we do not see any direct correlation between 
Peye and Prot when setting Prot = Pw These two parameters show 
random scatter with a Pearson correlation coefficient r = -0.27. 
Thus, we decided to search for a connection of Piot/Pcyc to Ro”' 
and logP^jjj. 

We define the inverse Rossby number, equivalent to the 
Coriolis number, as Ro ~* = Co = 2Q.Tc - AnTdPmi 
(iBrandenburg et al]ll998h . where Tc is the convective turnover 
time. The convective turnover time cannot be observed directly 
and needs to be determined in some other way. We used the 


Article number, page 19 ofl24l 























































A&A proofs: manuscript no. solan 



log Ro ^ l°g R HK 

Fig. 14. Left: log P^ilPcyc vs. logRo^*. Our stars are marked black if their cycles are graded “fair” or better and white if their cycles are “poor”. 
Triangles denote star s with cycles labelled “long”, while diamonds are used for the rest of the sample stars. Grey points denote reference data from 
ISaar & Brandenburg lll999h . Vertical dashed lines connect the long and short cycles found on the stars HD 63433, HD 82443, and HD 116956. 
The curved dashed line is a Gaussian process fit to the activity branches apart from the inactive branch, approximately located above the light grey 
dashed line. Approximate ranges of the active, transitional and superactive branches are denoted by the horizontal arrows labelled “A”, “T” and 
“S”. Right: Same as on the left but for log PmJPcyz vs. log RJjk- Stars in our sample that have detected active longitudes are circled in the plot. The 
grey vertical dashed line marks the approximate divide between the active stars with active longitudes and the less active ones without them. The 
scale for log P^aJPcyc is identical in both panels. 


semi-empirical formula of iNoves et al.l (1 19841) which gives Tc as 
a function of the B - V colour. The resulting Tc and logRo^* 
values are listed in Table |3 

When the logProt/Rcyc values are plotted against logRo^* 
(Fig . [141 left panel) our results align nicel y along the active and 
transit ional branches defined by the data of ISaar & Brandenburd 
(Il999l) . One of our stars (HD 73350) diverges from this pattern 
and lies closer to their inactive branch near the upper left corner 
of the plot. Our data shows that the superactive branch at high 
logRo ' and low logProt/Rcyc has the transitional branch as its 
continuation. This in turn meets smoothly with the active branch 
at low log Ro ' and high log Prot/Rcyc forming a long and mean¬ 
dering logPiot/Rcyc sequence. The approximate logRo ' ranges 
of the active, transitional and superactive branches on this se¬ 
quence are denoted by the atTows labelled “A”, “T” and “S” at 
the bottom of the left panel in Fig. [14] 

Our stars occupy the region around the turning point between 
the active and transitional branches. This provides new data for 
locating the logRo * value at which the apparent change in the 
behaviour of the activity cycle lengths occurs. To determine this 
value we cho se our cycles graded “fair” or better and the cy¬ 
cle data from ISaar & Brandenbu^ d 19991) excluding the inac¬ 
tive branch. As a working criterion we defined all cycles with 
logProt/Rcyc > 0.5 log Ro ' - 2.6 as belonging to the inactive 
branch. This limit is marked by the light grey dashed line in 
the left panel of Fig[T4] We model this cycle data with a simple 
Gaussian process applying a Gaussian covariance function ove r 
log Ro * with a standard deviation of ctcov — 0.95 (iHara 3201 Ih . 
This produces a predictive mean for log Prot /Pcyc as a function of 
log Ro *, plotted as the black dashed curve in Fig. [14] The curve 
peaks at logRo * = 1.42, locating the turnoff point between the 


ac tive and the transitional bran ches near this value. As discussed 
bv ISaar & Brandenburd (1 19991) . this point may mark a transition 
from antiquenching of the a-effect on low values of log Ro * and 
magnetic field strength to the more traditional quenching on the 
higher values. 

One problem with comparing Piot/Pcyc to Ro ' is that the 
Rossby number cannot be observed directly and is itself a func¬ 
tion of Prot- Some level of spurious correlation between the two 
values is thus expected. A resolution t o this problem is t o use 
the relation between Ro ' and logP[jj^ (iNoves et al.lll984l) . The 
chromospheric emission index logP[jj^ can be observed inde¬ 
pendently of Pfot and Pcyc and be used to estimate the spuri¬ 
ous correlation. The right panel of Fig. [14] shows our results 
plotted as log Pmf /Pcyc vs. l og over the reference data from 
ISaar & Brandenburd (1 1999h . A logP[jj^ value is not available for 
all of the stars in the reference data so the plot appears sparser 
than the logProt/Pcyc vs. logRo^’ plot. Nevertheless, the same 
active and inactive branches and the turnoff to the falling tran¬ 
sitional branch can be recognized f rom it. Now the Vaughan- 
Preston gap at lo gP(„. » -4.75 (IVaughan & Preston] Il980t 
iHenrv et^l 19961) approximately divides the active and inactive 
branches from each other. The similarity in the structures seen in 
the two plots gives a reassuring picture that the activity branches 
seen in the logProt/Pcyc vs. logRo ' plot are not dominated by 
the suspected spurious correlation. 

A feature that is n ot seen in the cycle data of 
ISaar & Brandenburd (1 19991) but is evident in our data, especially 
in the logProt/Pcyc vs. logPJ^j^ plot, is the split of the activ¬ 
ity cycles into two parallel sub-branches both within the active 
and the transitional branches. The sub-branches are in fact de- 
hned so well in our data that even the uncertain “poor” cycles 
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follow them closely. We find that none of our stars falls be¬ 
tween the sub-branches. It is also significant that for the three 
stars with double cycles the two cycles fall on the separate sub¬ 
branches, the shorter on the upper and the longer on the lower 
sub-branch. Towards higher logRo^' or logR^jj, the behaviour 
of the sub-branches is less certain. However, if the “poor” cycle 
of HD 29697 at high log can be seen to fall on the upper sub¬ 
branch, this sub-branch continues as far down the transitional 
branch as the lower sub-branch does. 

It is unclear how the sub-branc hes should be r e lated to the 
previous results by other authors. iBohm-Vitensel (l2007h sug¬ 
gested a division of the active branch into two sub-branches but 
these cannot be easily interpreted as our sub-branches. 

7. Conclusions 

In this paper we have analysed between 16 and 27 years of dif¬ 
ferential photometry from 21 young solar-type stars and char¬ 
acterized the nature of their activity with the aid of time series 
analysis. Furthermore, we have obtained high resolution spectra 
of the stars and calculated new chromospheric emission indices 
for them. The results show that the activity related phenomena 
of these stars follow a number of different trends. 

We estimated the differential rotation of the stars by de¬ 
termining the range of the observed photometric period vari¬ 
ations following the methodology of Lletsul (Il993l) . There are 
certain problems in interpreting the Z parameter value directly 
as a measure of the differential rotation since both active re¬ 
gion growth and decay and computational instabilities can in¬ 
crease the uncertainties in the measured period values. We con¬ 
clude, however, that these effects do not dominate our results. 
The remaining uncertainty concerns the unknown latitude ex¬ 
tent of the spot activity on the observed stars and how the pe¬ 
riod variation ranges should be scaled into the differential rota¬ 
tion coefficient k. Nevertheless, the measured period variations 
correlate with the surface differential rotation and reveal on the 
one hand a steep trend of increasing k towards longer P^t, as 
k oc PjQj®, and on the other a much flatter dependence of AQ 
on Q, as AQ oc q-o 36 pf these results agree qualitatively 


with previous studies (iHenrv et al. 

1995al:lDonahue et al.lll996l: 

iBarnes et al.ll2005t iReinhold et al. 

20131). However, it has to be 


noted that the published results for the surface shear are widely 
scattered around constant AQ, and it is not immediately clear 
how its underlying rotation dependence should be interpreted. 
We did not find any temperature dependence of the differential 
rotatio n, w hich again fits the rec ent results of iKiiker & Riidig^ 
(120111) and iReinhold et al.1(l2Q13l) for stars with effective temper¬ 
atures between 3500 K and 6000 K. 

It should be borne in mind that there are still considerable 
uncertainties in measuring differential rotation with the Z pa¬ 
rameter. It is possible to obtain information on the rotation and 
temperature dependence of k and AQ using Z for a population 
of stars, but the differential rotation estimates of individual stars 
can have substantial errors and should be taken with a grain of 
salt. 

Activity cycles were found from the photometry of nearly all 
of our stars with some level of significance (see Table [3]l. When 
the ratio of the rotation and cycle periods was plotted against 
log Ro * or logR^jj,, our stars fell n eatly on the active and transi- 
tional activity branches defined by ISaar & Brandenburg (1 19991) 
with one star being closer to their inactive branch. Our stars 
show that the transitional branch meets smoothly with the ac¬ 
tive branch at around logRo ' = 1.42. Thus, the superactive 
and transitional branches form a continuous trend of decreasing 


scaled cycle periods with decreasing Ro * down to the turnoff 
point of meeting with the active branch and turning smoothly 
into an opposite trend. Something not clearly seen until now is 
the split of the active and transitional branches into two parallel 
sub-branches, as is visible in Fig. [14] This finding suggests that 
there are multiple cycle modes available for the dynamos across 
a wide range of rotation rates. 

Active longitudes were also frequently found, although they 
do not appear to be quite as common as the activity cycles (see 
Table |3]l. We found temporary or persistent active longitudes 
from 11 of our stars, all of them belonging to the more active 
part of our sample. There appears to be a divide at approximately 
logRJ^j^ = -4.46 so that nearly all of the stars on the less ac¬ 
tive side of the limit have no active longitudes, while nearly all 
on the more active side have them. This suggests that the large 
scale magnetic fields are dominated in the less active stars by 
axisymmetric dynamo modes up to a certain limit, after which 
the more active stars develop strong non-axisymmetric modes. 
In Fig. [T4| the limit appears close to the break point between 
the active and transitional activity branches, suggesting a deeper 
connection between the opposite trends seen in the activity cycle 
lengths and the dominant magnetic field geometry. 

In all but one case the estimated active longitude periods Pal 
were shorter than the mean photometric rotation periods P^. 
In this single deviating case of opposite order of Pai and 
the two estimated periods were practically identical. We com¬ 
pared the two periods on each star and found that for seven 
of them a single underlying rotation period does not seem ad¬ 
equate to simultaneously describe them. Our interpretation is 
either that the magnetic structures supporting the active longi¬ 
tudes in these stars reside on deeper levels in the interiors, where 
the rotation is faster or that these stars have prograde propagat¬ 
ing azimuthal dynamo waves resulting in a trend of systemati¬ 
cally shortened active longitude periods. Similar migrating ac¬ 
tive l ongitudes have been observ ed previously from active giant 
stars (iHackman et al.l2()Tlll2013l) and propagating azimuthal dy¬ 
namo wayesj^ebeen seen in direct numerical dynamo simula¬ 
tions (ICole et al.ll2014h . The picture is somewhat complicated by 
the observation that the active longitudes commonly show some 
level of migration around the average active longitude period. 
Both the migration patterns and the disappearance and reappear¬ 
ance of the shorter-lived active longitudes show little regularity 
and still need to be fully explained. 
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Table 5. Rotation periods, activity indices, age estimates and identifications of kinematic group membership for the sample stars from previous studies. 
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log/?' 

log/?x 

age“ [Myr] 

group® 

Star 

Pmt[d] 

log/?' 

log/?x 
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group® 

HD 1405 

1.745^ 

-3.85^® 


20® 

49,50 
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1.76159^ 
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Table 5. continued. 


Star 

Pmt[d] 

log/?' 
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age“ [Myr] 

group® 
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Notes. Age estimates are based on: Tiso isochrone fitting, Tu lithium abundance, Tgyro gyrochronology, Thk Can H&K emission, Tx X-ray emission. ® Labels for kinematic groups: "AB" AB 
Dor Moving Group, "HLA" Her-Lyr Association, "Hya" Hyades Supercluster, "IC" IC 2391 Supercluster, "LA" Local Association, "Pie" Pleiades Moving Group, "UMa" UMa Moving Group, 
"THA" Tuc-Hor Association, "TWA" TW Hya Association 
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